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6.1 


ABSTRACT 

Fluctuating  temperature  structures  in  the  atmosphere 
induce  phase  perturbations  in  a  propagating  laser  beam.  These 
turbulent  conditions  occur  throughout  the  atmosphere  and  cause 
the  laser  beam  to  spread  and  alter  its  centroid.  There  are 
several  methods  to  measure  the  parameters  of  optical 
turbulence  in  the  atmosphere,  but  few  that  will  determine  them 
as  a  function  of  altitude  at  all  levels.  One  method  of 
measuring  altitude  profiles  of  turbulence  is  with  a 
temperature  probe  launched  via  a  balloon  system. 

This  thesis  involves  the  development  of  a  differential 
temperature  probe  sensor  to  measure  the  temperature 
fluctuations  at  all  altitudes  in  the  atmosphere.  In  addition, 
it  investigates  the  effect  of  solar  heating  on  the  probes  in 
the  atmosphere  and  the  subsequent  effects  on  the  measurements. 
A  validation  of  the  probe  system  was  made  by  a  comparison  test 
with  an  acoustic  echosounder  developed  earlier.  In  addition 
to  validating  the  probe  system,  the  absolute  CT2  analysis  of 
the  echosounder  was  confirmed. 
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I .  INTRODUCTION 

Atmospheric  conditions  will  cause  severe  degradation 
along  the  optical  path  of  a  ground  to  space  weapon  or 
surveillance  laser.  There  are  several  causes  for  this 
degradation,  they  are  1)  absorption  and  scattering  by  aerosols 
such  as  rain  and  clouds,  2)  distortion  by  thermal  blooming 
and  3)  distortion  by  atmospheric  turbulence  [Ref.  1]. 
Absorption  and  scattering  can  be  controlled  by  varying  the 
wavelengths  of  the  laser  and  having  multiple  sites  to  insure 
at  least  one  has  a  cloud  free  line  of  sight.  Thermal  blooming 
is  the  heating  of  the  medium,  through  which  the  laser  beam 
propagates.  It  is  due  to  absorption  of  the  radiation  by 
molecules  and  aerosols  and  the  consequent  distortion  of  the 
beam  due  to  density  reductions  brought  on  by  the  heating. 
Choosing  approximate  wavelengths  of  the  laser  which  have  low 
atmospheric  absorption  in  the  atmosphere  reduces  this  effect. 
Atmospheric  turbulence  is  difficult  to  deal  with  since  there 
is  no  way  to  avoid  it. 

The  major  effect  of  turbulence  on  an  optical  beam  is  the 
limitation  of  the  mutual  coherence  lengths.  For  example,  an 
average  ground  to  space  coherence  length  for  the  atmosphere 
is  on  the  order  of  5  cm  therefore  a  ground  based  laser  having 


a  4  meter  diameter  mirror  will  deliver  less  then  1/1000  of  its 
original  power  onto  a  target.  Adaptive  optics  provides  a  means 
for  reducing  these  turbulent  effects.  It  corrects  for  the 
effects  of  turbulence  by  altering  the  wavefront 
characteristics  of  a  beam  [Ref.  2]  using  deformable  mirrors 
or  nonlinear  optical  materials  such  as  Barium  Titanate. 

Measurements  of  the  turbulence  from  the  surface  to  an 
altitude  where  the  turbulence  is  no  longer  significant 
(approximately  3  0  km.)  are  needed  for  several  reasons.  The 
most  important  is  to  be  able  to  characterize  the  turbulent 
profiles  of  the  atmosphere  at  different  locations  in  order  to 
select  the  best  site  for  a  ground  based  optical  or 
surveillance  system.  Each  site  will  have  a  different 
turbulence  profile  since  it  depends  on  the  local  geography  as 
well  as  the  upper  atmospheric  conditions  controlled  by  general 
meteorological  patterns.  Another  reason  for  the  measurements 
is  to  determine  vertical  distribution  of  the  turbulence 
parameters,  such  as  the  coherence  length,  which  affect  the 
design  of  adaptive  optics  systems,  or  the  signal  processing 
transformations  in  a  surveillance  system. 

There  are  several  instruments  for  measuring  turbulence 
parameters,  some  of  which  measure  the  index  of  refraction 
structure  parameter  and  others  that  measure  the  temperature 
structure  parameter.  Some  of  the  methods  include  analysis  of 
star   trails   on   photographic   emulsions   [Ref.   3],   an 


isoplanometer  which  measures  the  isoplanatic  angle  through 
stellar  scintillation  [Ref.  4],  and  a  Modulation  Transfer 
Function  (MTF)  device  for  determining  the  coherence  length 
[Ref.  5].  These  devices  measure  important  integrated 
parameters  but  they  cannot  measure  the  vertical  profile  of 
turbulence.  An  acoustic  echosounder  which  is  similar  in  design 
and  construction  to  a  SONAR  [Ref.  6],  measures  a  vertical 
profile  but  is  usually  limited  to  several  hundred  meters 
range.  Greater  vertical  resolution  occurs  at  the  expense  of 
decreased  maximum  range. 

In  order  to  get  a  measurement  of  the  vertical  profile  of 
turbulence  to  20-50  km  a  device  must  be  able  to  be  launched 
on  a  sounding  balloon  or  an  aircraft.  An  example  of  this  type 
of  device  is  the  thermosonde  originally  designed  by  GTE 
Sylvania  and  revised  by  the  Air  Force  Weapons  Laboratory  and 
Tri-Con  Associates  Inc.  and  built  and  used  by  the  Air  Force 
Geophysics  Laboratory  [Ref.  7]. 

This  thesis  is  an  attempt  to  design  and  build  a 
temperature  sensing  probe  system  to  measure  the  vertical 
profile  of  the  temperature  structure  parameter  which  is  a 
measure  of  the  vertical  profile  of  turbulence.  Although  such 
systems  exist,  they  are  1)  expensive,  greater  than  $2  000  per 
launch  and  2)  require  extensive  calibration.  The  purpose  of 
this  thesis  was  to  develop  an  inexpensive  device  (the  system 
cost  is  approximately  $150  per  launch)  which  is  simple  to 


operate  (it  is  a  self-calibrating  device) .  Additionally  the 
effects  of  solar  heating  of  the  probes  are  investigated. 

The  results  of  a  comparison  with  a  well  developed 
acoustic  echosounder  indicates  the  system  will  be  effective 
in  measuring  turbulence  in  the  atmosphere.  The  studies  of  the 
solar  heating  effects  indicate  the  only  source  of  error  due 
to  solar  heating  would  be  from  the  sun/shade  effect  on  the  two 
probes  and  if  that  is  corrected  for,  it  will  be  accurate  up 
to  3  0  km. 
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II.   BACKGROUND 

A.   TURBULENCE  PARAMETERS 

Small  temperature  variations  carried  by  the  turbulent 
velocity  field  in  the  atmosphere  produce  small  phase 
perturbations  in  an  optical  plane  wave  propagating  through  it. 
These  perturbations  randomly  distort  and  convolve  the  phase 
of  a  plane  wave.  There  are  three  atmospheric  parameters  which 
must  be  determined  prior  to  any  attempt  made  at  compensating 
for  these  atmospheric  distortions.  These  parameters  are  the 
refractive  turbulence  structure  parameter,  CN  ,  the  spatial 
coherence  length  of  the  atmosphere,  r0  ,  and  the  isoplanatic 
angle  9 0. 

The  most  important  of  these  parameters  is  CN2.  Tatarski 
[Ref.  3]  states  that  one  way  to  deal  with  a  non-stationary 
problem,  which  includes  all  atmospheric  parameters,  is  to 
define  a  function  in  terms  of  a  difference,  he  than  defines 
the  structure  function  for  index  of  refraction  as, 

Dn(iV?2)  =<  [N(?2)-N[?1)]2>  ,  (i) 

where  <  >  denotes  an  ensemble  average.  If  we  assume  the 
atmosphere  to  be  homogeneous  and  isotropic  over  small  regions 
the  structure  function  can  be  rewritten  as, 


DN(r)  =  <[N(r2)-N(r1)]2>/  (2) 


where  r  is  r,,-r2.  By  dimensional  analysis  Kolomogorov  [Ref.  3] 

•  2/3 

showed  that  the  structure  function  has  an  r    dependency. 

•  •  2    •        2/3 

Consequently  DN  is  proportional  to  a  constant  CN  times  r  , 


DN  =  CN2r2/3.  (3) 


CN  is  the  refractive  turbulence  structure  parameter,  a  mean- 
square  statistical  average  of  the  difference  in  the  index  of 
refraction  between  two  points  separated  by  r12, 


CN2  =  <(N2-N1)2>/r122/3.  (4) 


2/3  •      •  • 

The  r  normalization  extends  from  an  inner  scale  10,  on  the 
order  of  millimeters,  to  an  outer  scale  of  L0,  on  the  order  of 
meters  [Ref.  3].  These  fluctuations  in  the  index  of  refraction 
arise  from  variations  in  density  caused  by  temperature 
fluctuations  in  the  turbulent  velocity  field. 

These  density  variations  in  the  atmosphere  alter  the 
phase  of  an  optical  beam  being  propagated  through  it.  The 
Optical  Transfer  Function  (OTF)  characterizes  the  integrated 
phase  perturbations  of  an  optical  beam.  It  is  a  measure  of  the 
correlation  of  the  electric  fields  of  the  optical  beam 


perpendicular  to  the  direction  of  propagation.  Although  the 
atmosphere  is  not  homogeneous  or  isotropic,  Tatarski  [Ref.  8] 
postulates  the  idea  of  local  homogeneity  and  isotropy,  which 
states  that  over  some  region  R,  comparable  to  the  outer  scale 
length  LQ/  we  can  assume  the  atmospheric  random  variables  are 
homogeneous  and  isotropic.  The  modulus  of  the  OTF  is  the 
atmospheric  Modulation  Transfer  Function  (MTF) .  Fried  [Ref. 
9]  introduces  the  parameter  r0  to  characterize  the  MTF.  It  is 
related  to  the  refractive  turbulence  structure  parameter  by. 


ro  =  2-1 


1.46k2 


rL  0    -,-3/5 

C  2(z)dz 
0  n 


(5) 


where  r0  is  the  spatial  coherence  length,  k  is  the  wave 
number  (2*/\) ,  and  CN2  is  the  refractive  turbulence  structure 
parameter  along  the  optical  path  of  length  L  [Ref.  5] . 

The  other  measure  of  spatial  coherence  in  the  atmosphere 
is  the  isoplanatic  angle  6 0.  It  is  similar  to  r0  in  that  it 
is  the  dependence  of  the  optical  transfer  function  of  a  system 
for  different  angles  to  the  source.  The  two  parameters  r0  and 
60  are  conjugate  pairs.  $0  looking  up  is  equivalent  to  r0 
divided  by  the  path  length  looking  down  [Ref.  4]  and  vice 
versa.  A  more  formal  definition  is  that  $0  is  an  angular 
measure  of  spatial  coherence,  it  is  the  limiting  angle  for 
which  an  electromagnetic  wave  from  a  source  will  follow  the 
same  optical  path  length  to  a  receiver. 
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Figure  1.  Comparison  of  LASER  Beam  Propagating 

Through  Vacuum  and  Atmosphere,  showing  rQ, 
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Figure    2.    Two   Sources   Showing   the   Effects    of 
the   Isoplanatic  Angle. 


If  we  consider  two  paths  through  the  turbulence,  the 
isoplanatic  angle  relates  the  mutual  coherence  e'1  point  of  the 
E  field  between  the  two  paths.  Fried  [Ref.  9]  expresses  the 
isoplanatic  angle's  dependence  on  CM2  as, 


°o  = 


*L_    9.        S/3.    1~3/5 

0 


2.905k2  [  C  2(z)z5/3dz 
Jn    n 


(6) 


where  z  is  the  altitude.  Note  the  z5/3  spherical  weighting 
factor  that  emphasizes  turbulence  far  from  the  optical  system. 
The  preceding  paragraphs  clearly  show  the  importance  of 
the  CN2  parameter,  it  not  only  defines  the  measure  of 
turbulence,  it  also  determines  the  spatial  coherence  r0  and 
the  isoplanatic  angle  6  0.  The  problem  lies  in  that  high 
resolution  profiles  of  CM2  are  difficult  to  measure,  it 
requires  complex  detectors  and  optical  imaging  systems. 
Instead  we  can  define  a  temperature  structure  parameter  CT2 
similar  to  C 2  where, 


CT2  =  <(T2-T1)2>/r122/3,  (7) 


which  can  be  measured  by  several  different  methods.  The 
fluctuations  of  the  index  of  refraction  are  due  to 
fluctuations  in  the  density  of  the  atmosphere  and  if  we  can 
assume  that  the  density  fluctuations  are  due  solely  to 
temperature  fluctuations,  then  CT2  is  related  to  CN  by, 
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cn2  = 


79xlO"6P 


T2 


^   ,  (8) 


where  P  is  the  atmospheric  pressure  in  millibars  and  T  is  the 
atmospheric  temperature  in  Kelvins  [Ref.  10].  The  assumption 
that  index  of  refraction  fluctuations  are  due  only  to 
temperature  fluctuations  and  that  humidity  fluctuations  are 
insignificant,  is  valid  when  the  Bowen  Ratio  B  (ratio  of 
sensible  heat  flux  to  latent  heat  flux)  is  greater  than  0.3 
[Refs.  11,12].  Below  this  value,  humidity  fluctuations  are 
significant. 

B.   SCALE  LENGTHS 

In  Kolomogorov's  definition  of  the  structure  function  he 
assumed  local  homogeneity  over  a  region  bounded  by  the  inner 
and  outer  scale  lengths  10  and  L0.  These  scale  lengths  vary 
from  hundreds  of  meters  at  the  outer  scale  lengths  down  to 
millimeters  for  the  inner  scale  lengths  [Ref.  13].  The  outer 
scale  length  is  the  size  of  the  turbulent  fluctuations  at  the 
onset,  while  at  the  inner  scale  viscosity  dissipates  the 
energy  of  turbulence  as  heat.  Kolomogorov  called  the  region 
between  10  and  L0  the  inertial  subrange  and  as  long  as  the 
distance  r,  in  Equations  (4)  and  (7),  is  within  this  inertial 
subrange  these  equations  are  valid.  Kolomogorov  expressed  the 
power  spectral  density  of  the  turbulence  in  this  region  by, 
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<t>(K)  =0.033Cn2K-H/3  f  (9) 


where  2*L„   «  K  «  2*1/  .  Von  Karman  took  this  definition 

0  0 

further  to  include  the  ranges  for  eddy  sizes  greater  than  L0, 


r(  11/6)  tt-9/2   9   , 

♦n^  ~  r(l/3)   ~8—  <^n>Lo 


K2 

1+  K 


\2 


11/6 


(10) 


where  K0=2wL0"1  and  <5n2>  is  the  variance  of  the  refractivity 
fluctuations  and  is  related  to  CN2  by, 


CN2  «    1.9<5n2>K02/3,  (11) 


[Ref.  10].  Figure  3  shows  the  spectrum  of  turbulence  for  all 
scale  lengths,  the  inert ial  subrange  shows  Kolomogorov's 
linear  description  of  turbulence  while  in  the  regions  above 
2wl0*1  viscosity  effects  dominate  and  below  2*L0~1  Von  Karman 's 
spectrum  defines  the  turbulence  [Ref.  14]. 

It  is  important  to  understand  the  effects  of  the 
turbulence  spectrum  in  all  three  regions  because  of  the  errors 
introduced  due  to  the  lack  of  correlation  in  the  regions  above 
and  below  the  inertial  subrange.  Care  must  be  taken  in 
choosing  the  correct  value  of  r  to  cover  the  scale  lengths 
inside  the  inertial  subrange.  Additionally  the  temporal 
frequency  of  the  turbulent  fluctuations  is  related  to  the  wave 
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number  K  and  the  wind  speed  moving  the  turbulence  past  the 
probes.  Therefore  to  measure  the  thin  transition  layers  in  the 
atmosphere  accurately  [Ref.  15]  the  errors  due  to  various 
scale  sizes  must  be  determined.  These  calculations  will  be 
carried  out  in  Chapter  4. 
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Figure  3.  Spectrum  of  Turbulence  Over  All  Scale 

Lengths . 


C.   SYSTEM  REQUIREMENTS 

A  fast  response  time  temperature  probe  attached  to  a 
sounding  balloon  system  can  measure  the  vertical  profile  of 
the  temperature  structure  parameter.  The  design  of  the 
temperature  probe  system  provides  a  simple,  low  cost  method 
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that  is  capable  of  resolving  the  thin  stratified  layers  of  a 
stable  atmosphere  and  determining  the  temperature  gradient, 
thickness,  and  turbulence  of  these  layers.  In  designing  the 
system  the  key  points  of  consideration  were  1)  what  type  of 
probe  geometry  and  thereby  what  type  of  circuitry  and  2)  what 
type  of  temperature  sensing  element  to  use. 

The  probe  system  could  use  a  differential  measurement  or 
a  single  point  probe.  The  advantage  of  the  single  point  probe 
is  that  it  can  measure  CT2  from  either  the  variance  of  the 
data,  knowing  the  balloon  ascent  rate,  or  by  analyzing  the 
power  spectral  density. However,  this  would  require  a  high  data 
rate  (several  samples  per  second) ,  which  the  radiosonde 
systems  used  do  not  have.  Although  a  differential  system  would 
not  have  as  high  a  vertical  resolution  it  has  the  advantage 
of  providing  partially  reduced  data  which  relaxes  the  need  for 
a  high  data  rate.  Therefore  a  differential  system  greatly 
simplifies  the  telemetry  needed  for  the  system  at  the  expense 
of  more  complexity  in  the  sensor  itself .The  resolution  for  the 
system  would  still  be  satisfactory,  about  2  meters  of  vertical 
resolution  for  a  balloon  with  an  ascent  rate  of  2  meters  per 
second. 

There  are  several  different  choices  for  the  probes.  They 

i  can   be   made   from   resistance   wires,   thermistors,   or 

thermocouples.  The  resistance  wire  is  simply  a  fine  wire,  such 

as  platinum  or  tungsten,  with  a  known  resistance  as  a  function 
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of  temperature.  The  problem  with  a  resistance  wire  for  this 
application  is  that  it  requires  some  means  of  self-calibration 
due  to  its  dependence  on  Ohm's  law  and  temperature.  A 
thermistor  is  a  small  semiconductor  device  which  changes  its 
resistance  as  a  function  of  temperature.  It  has  a  larger 
change  in  resistance  vs.  temperature  then  a  metal  wire, 
although  it  is  non-linear.  However,  thermistors  are  large 
compared  to  a  fine  wire,  having  a  larger  thermal  mass  than  a 
probe  made  from  a  fine  wire,  which  increases  the  response  time 
and  the  susceptibility  to  solar  heating  of  the  device. 
Additionally,  both  the  thermistor  and  the  resistance  wire 
require  a  current  source  which  not  only  increases  the 
complexity  of  the  circuit  but  also  introduces  a  self  heating 
factor.  A  thermocouple  consists  of  two  fine  wires  made  of 
dissimilar  metals  welded  together,  which  produces  a  voltage 
difference  proportional  to  the  temperature.  For  an  in-depth 
discussion  of  thermocouples  and  the  thermoelectric  effect,  see 
Refs.  16  and  17.  Commercially  available  thermocouples  can  be 
made  of  very  fine  wire  (down  to  12.5  urn)  thereby  reducing  the 
thermal  mass  and  producing  a  faster  response  time.  A 
disadvantage  of  a  thermocouple  is  that  the  response  to 
temperature  is  small,  typically  40  ^V/°K.  This  places  severe 
requirements  on  low  noise  signal  processing.  An  operational 
requirement  for  this  type  of  device  would  be  the  knowledge  of 
the  mean  atmospheric  temperature  to  calculate  the  Seebeck 
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Coefficient,  which  is  the  derivative  of  the  thermal  emf  with 
respect  to  the  temperature.  Since  the  rawinsonde  system  used 
with  the  probe  provides  this  data,  it  is  easily  accomplished. 
Small  temperature  changes  of  1°  C  or  less,  are  expected  from 
the  two  probes  in  a  differential  system.  This  produces  a 
negligible  change  in  the  Seebeck  coefficient  therefore  no 
calibration  of  this  type  of  device  is  needed,  other  than 
knowing  the  gains  of  the  electronic  amplifiers. 
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III.  SYSTEM  DESIGN  AMD  DEVELOPMENT 


A.   AMPLIFIER  CIRCUIT 

The  probe  system  consists  of  a  pair  of  thermocouples 
connected  in  series  and  held  rigidly  1  meter  apart  by  an 
aluminum  tube.  An  amplifier  circuit  requires  high  voltage  gain 
and  very  low  noise  to  be  able  to  discriminate  the  signal  from 
the  background  due  to  the  low  voltage  produced  by 
thermocouples  (on  the  order  of  /xvolts  per  degree)  .  It  also 
requires  an  analog  root  mean  square  device  to  calculate  < (T2- 

Since  the  purpose  of  the  system  is  to  measure  temperature 
differences  on  the  order  of  hundredths  and  thousandths  of  a 
degree  and  the  thermocouple  output  voltage  is  about  40 
microvolts  per  degree  C  the  circuit  must  have  a  very  high  gain 
with  ultra  low  self  noise.  Therefore  the  circuit  must  be 
carefully  designed  and  built  to  reduce  any  sources  of  noise 
wherever  possible,  by  methods  such  as  matching  resistor  values 
as  exact  as  possible  and  placing  them  as  close  as  possible  to 
each  other  to  reduce  the  thermal  drift.  Other  examples  include 
thermal  insulation  for  the  entire  circuit  and  RF  shielding  for 
the  circuit  as  well  as  the  probe.  The  most  critical  component 
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of  the  circuit  is  the  ultra  low  noise,  precision,  high  speed 
op  amps  which  have  voltage  noises  less  than  those  of  50  ohm 
resistors. 

The  circuit  (designed  by  Prof.  Don  Walters  and  fabricated 
by  Dale  Galarowicz)  is  a  low  noise,  wide  bandwidth 
Instrumentation  amplifier.  Three  operational  amplifiers 
produce  a  gain  of  10,000  and  a  combination  low  and  high  pass 
filter  with  a  gain  of  5  producing  a  total  voltage  gain  of 
50,000.  The  circuit  contains  two  Linear  Technologies  LT1028 
ultra  low  noise  precision  high  speed  op  amps  in  the 
instrumentation  amplifier  portion  of  the  circuit.  These  op 
amps  have  a  gain  bandwidth  product  of  75  MHz  and  a  self  noise 
of  0.85  nV/Hz1/2  at  the  frequencies  desired.  Figure  4  shows  the 
noise  and  frequency  characteristics  of  this  op  amp.  The  filter 
for  the  circuit  uses  an  LT1057  op  amp  whose  noise  and 
frequency  response  characteristics  are  shown  in  Figure  5. 
Based  on  the  frequency  response  curves  in  Figures  4  and  5 
[Ref .  18]  and  the  high  pass  filter  the  circuit  has  a  frequency 
response  from  0.16  Hz  to  200  Hz.  The  other  major  component  of 
the  circuit  is  an  Analog  Devices  AD637  high  precision,  wide 
band  RMS-DC  converter.  It  has  a  bandwidth  of  600KHz  at  lOOmV 
RMS  and  an  averaging  time  constant  of  2  5  msec/micro  F.  The 
entire  circuit  is  powered  by  two  9  volt  dry  cells.  Figure  6 
is  a  complete  schematic  diagram  of  the  circuit. 
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Voltage  Noise  vs  Frequency 
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Figure  4.  Noise  and  Gain  Characteristics  of  the 
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The  circuit  has  to  be  well  shielded  against  RF 
interference  due  to  the  high  gain  of  the  circuit  and  the  fact 
it  will  be  operated  in  the  atmosphere  where  it  is  highly 
susceptible  to  all  types  of  RF  signals.  Additionally  the 
thermocouple  wires  must  be  stretched  out  over  a  distance  of 
1  meter  and  will  act  as  an  antenna.  The  aluminum  tube  shields 
them  from  RFI  except  at  the  two  endpoints,  greatly  reducing 
extraneous  signals.  The  circuit  is  inclosed  in  a  styrofoam 
casing,  to  reduce  thermal  gradients  across  the  circuit,  which 
is  then  covered  in  aluminum  foil.  The  circuit  itself  has 
ferrite  beads  and  an  LC  filter,  at  the  inputs,  to  further 
shield  the  op  amps  from  RFI  and  the  entire  circuit  is  built 
on  a  ground  plane  which  has  been  grounded  with  the  foil 
shield. 

The  result  of  the  aluminum  foil  and  styrofoam  shielding 
and  use  of  ultra  low  noise  op  amps  is  an  amplifier  capable  of 
measuring  the  extremely  small  voltages  produced  by  the 
thermocouple  probes.  The  circuit  was  tested  in  the  anechoic 
chamber  in  the  basement  of  Spanagel  Hall  to  minimize  any 
temperature  fluctuations  and  then  measurements  of  the  circuits 
self  noise  were  taken.  Figure  7  shows  the  noise  spectrum 
measured  by  a  Hewlett  Packard  HP3561A  Dynamic  Signal  Analyzer. 
The  large  noise  spike  below  1  Hz  is  due  to  1/f  or  flicker 
noise  which  occurs  in  all  amplifiers,  due  in  large  part  to 
surface  leakage  of  transistors  [Ref.  19].  Although  the  noise 
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normally  occurs  at  frequencies  up  to  100Hz  it  has  been  reduced 
by  the  use  of  a  high  pass  filter.  The  noise  spike  at  60Hz  is 
due  to  AC  powered  equipment  operating  in  the  anechoic  chamber. 
This  noise  should  vanish  when  the  system  is  used  in  the  field, 
since  it  is  powered  by  DC  batteries,  as  long  as  care  is  taken 
to  insure  it  is  not  near  a  large  AC  power  source.  The  plot 
clearly  shows  the  self  noise  output  of  the  circuit  is  well 
below  100  microVolts/Hz1/2  in  the  frequencies  of  interest. 
Since  the  circuit  gain  is  50,000  and  the  Seebeck  coefficient 
is  on  the  order  of  40  microvolts  per  degree  this  translates 
into  a  noise  induced  measurement  of  less  than  0. 00005°C/Hz1/2. 
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The  entire  system  was  set-up  and  run  in  the  anechoic 
chamber  with  caps  on  the  exposed  thermocouples  to  further 
reduce  the  signal.  After  running  for  about  2  hours  to  let  all 
air  currents  in  the  chamber  settle,  the  system  produced  the 
results  shown  in  Figure  8.  This  shows  the  noise  produced  by 
the  circuit  introduced  an  error  corresponding  to  a  CT2  of  less 
than  10"6,  which  is  two  orders  of  magnitude  lower  than  the 
lowest  CT2  needed  for  a  usable  probe  system.  These  results 
indicate  the  circuit  self  noise  is  well  below  that  which  would 
have  a  degrading  effect  on  the  results. 

The  entire  circuit  package  measures  3"  X  3"  X  3"  and 
weighs  less  than  six  ounces. 
B .   THERMOCOUPLE 

In  1821  Thomas  Seebeck  discovered  that  two  wires  of 
dissimilar  metals  joined  at  one  end  and  heated  produce  a 
voltage  difference  across  the  open  ends.  This  voltage  is  a 
function  of  the  junction  temperature  and  the  composition  of 
the  two  metals.  Since  that  time,  many  different  thermocouple 
types  have  been  produced  based  on  the  combination  of  the  two 
different  metals  used  and  having  different  thermoelectric  and 
physical  properties.  The  requirements  for  the  system  included 
use  in  the  atmosphere  from  the  surface  to  20  km,  therefore  the 
temperatures  vary  from  3  0°  C  to  -3  0°  C  with  as  large  a  Seebeck 
coefficient  as  possible.  The  Seebeck  coefficient  is  the  slope 
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of  the  voltage  versus  temperature  curve  at  a  given 
temperature.  Based  on  these  requirements  the  T  type  or  Copper- 
Constantan  thermocouple  was  selected.  An  alternative  choice 
would  be  the  E  type  or  Chromel-Constantan  thermocouple. 

The  Copper-Constantan  thermocouple  is  composed  of  a 
copper  wire  and  a  55%  copper  45%  nickel  wire.  It  has  a 
temperature  range  of  -200°  C  to  350°  C  and  is  suitable  for 
applications  where  moisture  is  present.  Table  1  shows  the 
thermoelectric  voltages  referenced  to  0°  C  for  a  Copper- 
Constantan  thermocouple  based  on  the  National  Bureau  of 
Standards  reference  tables. 

Table  1  shows  that  the  Seebeck  coefficient,  which  is  the 
unit  difference  in  voltage  for  each  temperature  change,  is  not 
linear  over  the  entire  temperature  range.  To  determine  the 
Seebeck  coefficient,  dV/dt,  valid  over  the  entire  temperature 
range  desired,  the  data  over  the  temperature  range  desired 
from  Table  1  was  plotted  and  then  a  polynomial  regression  was 
performed  to  find  the  equation  of  the  curve.  The  derivative 
of  the  curve  was  taken  to  find  the  Seebeck  coefficient.  A 
fifth  order  polynomial  fits  the  data  from  -100°  C  to  3  0°C.  The 
Seebeck  Coefficient  for  a  Copper-Constantan  thermocouple  is, 

dV/dt=3  .  8707xl0"2+8  .  5348xl0"5t-3  .  3135xl0"7t2 

-2.77432xl0"9t3-l.253xl0"1V,  (12) 
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Table   1 

VOLTAGES    FOR   A  TYPE   T  THERMOCOUPLE 

toe0               I               2               J               »               9               t  7               •               9             10  oca  t 

tlH«Moei  CCI«IC     VOlt«Gt     IN    6BSOIUU    Mill     VOltS 

-170  -6.15(  .{71 

-760  -6.2)7  -6.13*  -6.239  -6.761  -t.169  -6.261  -1,1)1  -6.1)3  -6.135  -4.13*  -».!)(  -1*0 

-1)0  -6.KI  -».l(7  -*.191  -».19(  -6. Jo*  -».20*  -4.21*  -1.111  -*.226  -l.lll  -1.1)1  -It) 

-no  -6.10)  -«.n«  -6.177  -t.no  -t.nt  -t.itt  -t.iti  -t.ito  -».  i » '  -».n*  -i.m  -no 

-1)0  -6.001  -1.011  -1.011  -».0>*  -».0*9  -6.0)9  -t.Oil  -6.07*  -6.01'  -6.096  -6.10J  -210 

-110  -5. ((9  -5.901  -3.91*  -3.92*  -3.9!(  -5.950  -5.962  -5.91)  -5.965  -5.996  -6.001  -110 

-110  -5.75!  -5.767  -5.7(2  -5.793  -3. 109  -5.(23  -».(!»  -3.(50  -3.(4!  -3.(7*  -».••»  -110 

-100  -5.603  -3.»l*  -3.636  -1.630  -3.665  -5.660  -3.695  -3.710  -3.72*  -3.739  -3.753  -100 

-190  -5.»J9  -5.636  -5.67!  -3.6(9  -5.306  -3.311  -3.339  -3.333  -3.571  -5.3(7  -5.603  -190 

-1(0  -5.161  -5.179  -5.197  -5.313  -3.153  -3.!5I  -5.!69  -3.5(7  -9.606  -5.611  -5.6)9  -no 

-170  -5.069  -5.069  -5.109  -1,111  -3.167  -3.1*7  -3.1(6  -3.205  -5.11!  -5.262  -5.161  -HI 

-160  -».(65  -«.((6  -6.907  -6.9M  -6.96(  -6.969  -6.9(9  -5.010  -5.030  -5.050  -5.069  -160 

-150  -».6»(  -6.610  -6.69)  -6.715  -6.757  -*.73(  -6.7(0  -6.(01  -6.(23  -6.(66  -6.(63  -150 

-1*0  -6.619  -6.662  -6.666  -6.6(9  -6.511  -6.3)3  -6.531  -6.5(1  -6.*03  -6.616  -6.66(  -I6| 

-1)0  -6.117  -6.202  -6.226  -6.131  -6.173  -6.299  -6.)2!  -6.)67  -6.171  -6.195  -6.619  -llo 

-WO  -1.921  -3.969  -).976  -6.000  -6.026  -6.051  -6.077  -6.102  -6.127  -6.152  -6.177  -120 

-110  -1.636  -1.6(6  -1.711  -1.717  -1.766  -1.791  -1.611  -1.(66  -1.(70  -1.(97  -1.911  -IK 
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where  t  is  in  °C  and  V  is  in  millivolts.  The  temperature 
difference  measured  by  the  probe,  is  the  output  voltage 
divided  by  the  circuit  gain  and  divided  by  dV/dt  determined 


at  the  specific  temperature.   For 


1°  change   in  air 


temperature  the  average  change  in  dV/dt  is  less  than  0.2 
microvolts  per  degree  C  or  less  than  0.5%  change  in  dV/dt. 
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Figure  10.  Schematic  of  Thermocouple  Probe. 

The  probe  consists  of  two  0.00254  cm  diameter  Copper- 
Constantan  thermocouples  held  rigidly  1  meter  apart  by  an 
aluminum  tube.  This  wire  has  a  resistance  of  998.3  ohms  per 
double  meter  and  the  circuit  requires  a  low  input  impedance, 
therefore  it  is  necessary  to  cut  the  fine  thermocouple  wire 
short  and  solder  larger  wires  to  them  to  run  the  distance 
between  the  thermocouples  and  the  circuit  input.  Copper- 
Constantan  wires  0.0254  cm  in  diameter  were  used  which  have 
a  resistance  of  9.983  ohms  per  double  meter.  It  is  essential 
to  solder  the  copper  to  copper  and  constantan  to  constantan 
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to  insure  no  other  thermocouple  junctions  are  formed  which 
could  cancel  any  voltage  signal  generated.  The  resulting  probe 
yields  an  input  resistance  of  30  ohms. 


a 
2  00  r 
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Figure  11.  Response  Time  Versus  Wind  Speed  of  a 
.  00254  cm  Copper-Constantan  Thermocouple 


The  reason  for  the  fine  wire  was  to  reduce  the  thermal 
mass  which  decreases  the  response  time  of  the  thermocouple  to 
temperature  change.  The  response  time  is  defined  as  the  time 
required  to  reach  63.2%  of  an  instantaneous  temperature 
change.  The  OMEGA  Handbook  [Ref.  20]  gives  the  response  time 
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for  a  0.00254  cm  thermocouple  as  0.05  sec  in  still  air  and 
0.004  sec  in  18.3  m/sec  air.  The  response  time  r  can  be 
related  to  the  velocity  by  the  equation, 

1   =  A  +  BJ—       ,  (is) 

T 

where  A  and  B  are  constants  for  the  specific  thermocouple 
material  size  and  are  determined  by  the  response  times  given 
in  Ref.  20  (for  this  case  A=20  and  B=53.783  sec*1).  This 
results  in  a  response  time  verses  velocity  curve  shown  in 
Figure  11.  From  this  curve  the  response  time  r  §  5  m/sec  is 
.0071  sec.  The  interaction  of  the  frequency  response  and  the 
inner  scale  size  will  be  discussed  in  the  Results  section. 

C.  COMPUTER  AND  CODE 

The  data  acquisition  and  processing  portion  of  the  system 
consists  of  a  Hewlett  Packard  model  217  computer  with  a  20 
megabyte  hard  drive  and  2  megabytes  of  memory.  It  also 
contains  an  Infotek  BC203  Basic  compiler,  an  Infotek  AD200 
analog  to  digital  converter,  an  HP9133  floppy  disk  drive,  a 
monitor  and  printer.  The  computer  receives  data  from  the 
circuit  via  the  analog  to  digital  converter.  It  reduces  the 
data  and  displays  a  CT2  verses  time  (altitude)  profile. 

The  analog  to  digital  converter  is  the  primary  component 
for  the  data  acquisition,  however  it  introduces  errors.  The 
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A  to  D  converter  has  a  small  DC  voltage  offset.  This  offset 
is  small  so  it  does  not  significantly  affect  the  results  until 
CT2  reaches  values  on  the  order  of  10"3,  as  is  discussed  in  the 
results  section.  Each  A  to  D  converter  has  a  different  offset 
so  each  time  a  different  computer  is  used  the  offset  must  be 
measured.  In  addition  the  noise  of  the  electronic  circuit 
produces  a  DC  offset.  To  find  the  total  DC  offset,  the  entire 
system  was  set  up  and  run  in  the  anechoic  chamber  to  insure 
no  external  signals  are  introduced.  The  software  has  a  feature 
for  inputting  the  offset  and  another  to  display  the  voltages 
directly  from  the  probe.  To  find  the  offset,  initially  set  the 
offset  to  zero  and  then  after  the  system  has  run  for  about  one 
hour,  to  settle  all  the  air  currents  generated  by  the  set-up, 
read  the  voltages  by  pressing  the  PRINT_RAW  key,  this  displays 
the  actual  voltage  offset  then  simply  reboot  the  system  and 
input  the  offset  read  from  the  raw  voltages.  The  DC  offset  for 
the  system  used  for  the  experiment  was  -5  millivolts. 

Appendix  A  summarizes  the  features  of  the  software  and 
lists  the  code  for  controlling  and  producing  output  from  the 
system.  It  is  almost  fully  automatic.  Once  it  has  started 
running,  the  only  corrections  to  be  made  are  to  update  the 
ambient  air  temperature  any  time  there  is  a  change  of  two  or 
three  degrees,  so  the  Seebeck  coefficient  can  be  recalculated. 
There  is  a  hard  function  key  for  updating  the  air  temperature. 
A  further  refinement  to  the  system  would  be  to  have  the 
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rawinsonde  temperature  measurement  device  automatically  update 
the  program  with  the  new  air  temperature  at  a  given  interval 
to  make  the  system  fully  automatic. 

D.   SOLAR  HEATING 

It  is  important  to  understand  the  effects  of  thermal 
radiation  on  the  probes  since  they  will  be  ascending  through 
the  atmosphere  both  during  the  day  and  at  night.  During  the 
day  radiative  heating  from  direct  and  reflected  solar 
radiation  will  heat  the  probes  above  ambient  air  temperature 
while  at  night  the  probes  will  be  cooled  due  to  Planck 
radiation  from  the  probes  to  space.  There  are  several  ways  in 
which  the  heating  or  cooling  of  the  probes  can  introduce 
errors  into  a  differential  measurement  device.  A  hot  wire 
anemometer  effect  can  introduce  errors.  Another  is  the 
difference  caused  by  one  probe  being  in  the  sun  and  the  other 
being  in  the  shade  or  at  night  where  one  probe  has  a  direct 
line  of  sight  to  the  earth  and  the  other  is  blocked  as  by  a 
cloud.  In  a  hot  wire  anemometer  effect  the  probes  are  heated 
above  the  ambient  temperature  and  then  velocity  fluctuations 
across  the  probes  vary  the  heat  transfer  rates  away  from  the 
probes,  creating  a  false  temperature  difference  [Ref.  21].  In 
the  sun-shade  effect,  as  the  probe  ascends  through  the 
atmosphere  one  end  may  come  into  shade  either  from  a  cloud  or 
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the  balloon  shadow.  The  probe  may  rotate  as  it  ascends  giving 
one  thermocouple  a  different  aspect  to  the  sun  than  the  other, 
creating  the  same  sun/ shade  effect  to  a  lesser  extent.  In 
either  case  one  thermocouple  will  receive  less  direct  solar 
radiation  than  the  other  thereby  changing  the  heat  flow  on  one 
thermocouple  but  not  the  other,  again  introducing  a  false 
difference.  The  same  applies  at  night  if  one  thermocouple  is 
shaded  from  the  earth  by  a  cloud  it  will  not  radiate  thermal 
energy  at  the  same  rate  as  the  other  will. 

The  net  heat  flow  to  or  from  a  body  in  the  atmosphere  is 
described  by  the  heat  transfer  equation, 

qnet=<3s  +  <3a  +  qt  +  3c  +  <3k  "  <3r  >  (14) 

where  qs  =  portion  of  direct  solar  radiation  absorbed 

qa  =  portion  of  the  atmospheric  radiation  absorbed 
qt  =  portion  of  terrestrial  radiation  absorbed 
qr  =  thermal  radiation  emitted  by  the  wire 
qk  =  net  conduction  to  the  wire  from  the  atmosphere 
qc  =  net  convection  to  the  wire  from  the  atmosphere 
[Ref.  22].    The  temperature  difference  can  be  determined  by 
first  setting  qnet  equal  to  0.  The  equation  can  then  be  reduced 
into  the  heating  and  cooling  portions  by, 

ESD    +    ESR    +    ELE    =    ELW   +    Ec    ,  (15) 
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where  Eg,  =  heating  due  to  direct  solar  radiation 

ESR  =  heating  due  to  solar  radiative  reflection  of 

the  atmosphere 
ELE  =  heating  due  to  long  wave  radiation  from  earth 
ELU  =  cooling  due  to  long  wave  radiation  from  wire 
Ec  =  convective  cooling 
[Ref.  23].  If  we  further  assume  the  thermocouple  to  be  a 
horizontally  oriented,  infinite  cylinder  with  the  top  half 
radiating  to  the  sky  and  the  bottom  half  radiating  to  the 
earth's  surface,   the  temperature  difference  between  the 
ambient  air  and  the  thermocouple  is  given  by, 


AT=   s 


n  a 


Rs  +  *<L 


a      g       4 


(16) 


where  h  =  average  convective  conductance 

cs  =  short  wave  emissivity  of  the  thermocouple 
cL  =  long  wave  emissivity  of  the  thermocouple 
a   =  albedo 

R$  =  short  wave  incoming  radiation 
R  =  long  wave  radiation  from  the  earth's  surface 
Ra  =  long  wave  atmospheric  radiation 
a      =  Stephan-Boltzmann  constant 
T  =  air  temperature 
where  h  is  defined  by  Kreith  [Ref.  22]  as, 
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h=Kk 

w 

where  D  =  wire  diameter 


n 
VD 


(17) 


V  =  wind  speed 

k  =  heat  conductivity  of  the  air 

v   =  kinematic  viscosity  of  the  air 
K&n  =  empirically  determined  dimensionless  constants 
based  on  the  Reynolds  number 
[Ref.  24]. 

The  Air  Force  Geophysics  Laboratory  [Ref.  23]  defines  the 
calculation  for  delta  T  in  a  similar  manner  based  on  the  same 
assumptions.  However,  when  comparing  the  two  forms  using  the 
same  parameters  (see  Figure  12)  there  is  a  considerable 
difference.  Review  of  both  treatments  shows  the  principle 
difference  lies  in  each  definition  of  the  convective 
conductance  h.  Campbell  [Ref.  24]  uses  Krieth's  form  of  h  (EQN 
17)  while  Brown  [Ref.  23]  assumes  an  average  value.  Another 
method  of  determining  h  is  with  the  Nusselt  number,  which  is 
a  dimensionless  number  used  in  describing  heat  transfer  and 
fluid  flows.  Kramers  [Ref.  25]  performed  extensive 
measurements  of  heat  transfer  on  spheres  and  cylinders,  from 
this  he  determined  the  Nusselt  number  to  be  a  function  of  the 
Reynolds  number  and  the  Prandtl  number,  another  dimensionless 
number  where, 
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Re  =  VDp/n    ,    and   Pr  =  c  u/k  ,  (18) 

where  D  =  diameter  of  the  cylinder 

V  =  velocity  of  the  fluid 

p    =   density  of  the  fluid 

/i  =  dynamic  viscosity  of  the  fluid 

k  =  thermal  conductivity  of  the  fluid 

c  =  specific  heat  of  the  fluid. 
From  all  the  available  data  he  showed  that  the  Nusselt  number 
for  a  cylinder  could  be  represented  by, 


Nu  =  0.9l(Pr)°-31(Re)0-385   , 0.1<Re<50  , 

and  (19) 

Nu  =  0.60(Pr)°-31(Re)0-50    ,  50<Re<10,  000. 


Based  on  the  Nusselt  number  the  convective  conductance  is, 

h  =  Nu*k/D.  (20) 

Campbell  included  some  experimental  results  in  his  paper. 
When  the  experimental  results  are  compared  with  Kreith's  and 
Kramers'  treatment  of  h  we  can  see  that  Kramers'  expression 
exactly  models  the  actual  data  (Figure  13) .  Using  Equation  20 
for  h  in  Equation  16  corrects  the  differences  seen  in  Figure 
12,  therefore  the  lower  curve  in  Figure  12  is  the  correct 
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model  for  solar  radiative  heating  in  the  atmosphere,  based  on 
the  results  obtained  using  the  experimental  data  of  Reference 
24. 
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Figure  12.  Comparison  of  AFGL  Model  and  Campbell's 
Model  of  Solar  Heating  of  Thermocouple 
Wires. 


A  discrepancy  noted  in  Campbell's  calculations  was  the 
values  used  for  short (visible)  and  long(IR)  wave  emissivities. 
Table  2  gives  the  emissivities  used  and  the  actual 
emissivities  from  a  1986  edition  of  the  CRC  handbook.  The 
corrected  values  were  used  to  recalculate  the  solar  heating 
and  the  updated  results  are  shown  in  Figure  14 . 
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Figure  13.  Comparison  of  Convective  Conductance. 


Table  2 

EMISSIVITIES  OF  CU-CONSTANTAN 

AND  TUNGSTEN 


Visible 

IR 

Campbell's 
Cu-Constantan 

.25 

.5 

CRC  Handbook 
Cu-Constantan 

.2 

.03 

CRC  Handbook 
Tungsten 

.5 

.03 
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Figure  14.  Corrected  Comparison  of  solar  heating. 

Based  on  Figure  14  it  is  now  possible  to  make  a 
determination  of  the  errors  introduced  by  solar  heating.  The 
structure  function  for  velocity  fluctuations  over  small  scale 
lengths  is, 


D(r)  =  3.83(er) 


2/3 


(21) 


where  c  is  the  dissipation  rate  [Ref.  26].  Actual  data  shown 
in  Reference  2  6  from  areas  of  highly  turbulent  velocity  fields 
indicates  the  average  dissipation  rate  is  on  the  order  of 
3xl0"5  m2  sec"3.  For  the  probe  system  with  r  equal  to  1  meter 
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this  yields  velocity  fluctuations  on  the  order  of  0.06  m/sec. 
Figure  15  shows  the  solar  heating  errors  due  to  the  hot  wire 
anemometer  effect  with  velocity  fluctuations  of  this 
magnitude.  Figure  16  shows  the  differences  between  the  curves 
of  Figure  15.  It  indicates  the  temperature  differences  are 
negligible  (on  the  order  of  .001  degrees  or  less). 

The  shading  effect  can  be  determined  by  eliminating  the 
direct  solar  radiation  component  in  the  equation.  This  will 
show  the  maximum  error,  if  it  is  a  matter  of  the  probe 
changing  aspects  to  the  sun  the  errors  will  be  proportionally 
less.  If  direct  solar  radiation  is  completely  removed  there 
will  be  very  slight  heating  of  the  probe,  due  to  incoming 
terrestrial  radiation.  The  temperature  difference  between  the 
two  probes  will  be  approximately  equal  to  the  amount  of  solar 
heating  on  one  probe  as  seen  in  Figure  14 .  This  indicates  a 
major  source  of  error  since  the  temperature  difference  at 
higher  altitudes  is  on  the  order  of  0.2  °C  and  the  balloon  can 
rotate  as  it  ascends. 

Based  on  these  calculations  if  some  method  is  devised  to 
limit  the  rotation  of  the  probe  to  eliminate  the  sun/ shade 
effect,  the  probe  system  can  measure  CT2  values  accurately  up 
to  3  0  km  altitude  without  significant  errors. 

Reference  7  described  a  CT  thermosonde  used  by  the  Air 
Force  Geophysics  Laboratory.  Data  measured  by  this  system 
indicates  an  order  of  magnitude  jump  in  the  values  of  CN2 
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just  after  sunrise  [Ref.  27].  This  increase  appears  to  be  an 
artifact  of  the  instrument  rather  than  actual  turbulent 
processes  due  to  the  fact  it  occurs  so  rapidly.  Since  it 
occurs  at  sunrise  a  logical  assumption  is  that  it  is  due  to 
solar  heating,  therefore  a  great  deal  of  time  has  been  spent 
in  determining  these  effects.  Figure  17  shows  the  hot  wire 
anemometer  effect  on  a  4  micron  tungsten  resistance  wire.  The 
values  for  the  emissivity  of  Tungsten  are  taken  from  Table  2. 
This  indicates  solar  heating  does  not  effect  the  measurements 
of  CT2  consequently  the  rise  in  the  value  of  Cn2  at  sunrise  may 
be  actual. 
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Figure  17.  Hot  Wire  Anemometer  Effect  on  4  micron 

Tungsten  Wire. 


40 


IV.  RESULTS 

A.   EXPERIMENTAL  PROCEDURE 

The  experimental  measurements  served  two  purposes,  first 
they  were  carried  out  to  validate  the  probe  measurement  system 
and  second,  they  were  used  to  validate  the  CT2  measurement 
capabilities  of  the  acoustic  echosounder [Ref s.  6,28].  The 
acoustic  echosounder  calculates  a  15  minute  time  averaged 
value  for  the  temperature  structure  parameter  as  a  function 
of  altitude,  however  uncertainties  in  the  antenna  beam  shape, 
side  lobes,  transducer  efficiencies  and  atmospheric 
attenuation  produce  uncertainties  in  the  absolute  value  of  CT2 
calculated  by  the  echosounder  [Ref.  28].  Therefore  independent 
verification  of  the  CT2  values  must  be  made  to  validate  the 
acoustic  echosounder.  If  both  values  agree  this  is  a  positive 
indication  that  both  systems  are  measuring  accurately. 

The  acoustic  echosounder  is  a  high  frequency  device  which 
analyzes  the  atmospheric  density  fluctuations  within  the  first 
200  meters  of  the  atmosphere.  The  echosounder  consists  of  a 
Hewlett  Packard  HP  217  computer  to  control  the  system  and 
acquire  and  reduce  the  data,  an  HP  3314A  function  generator 
which  produces  the  pulsed  signals,  an  amplifier  and  the  array 
of  speakers  which  acts  as  a  transmitter/receiver.  The  system 
operates  at  5KHz  and  produces  a  100  cycle  sinusoidal  burst  of 


41 


18  acoustic  watts.  The  antenna  array  consists  of  19 
piezoelectric  speakers  in  a  close-packed  hexagonal  array 
enclosed  in  a  55  gallon  plastic  trash  container  lined  with 
lead  and  foam  to  reduce  side  lobes  as  much  as  possible.  The 
minimum  range  of  the  device  is  approximately  six  meters  based 
on  the  recovery  or  "ring"  time  of  the  speakers  and  the  maximum 
range  is  about  200  meters  based  on  the  frequency  used.  Figure 
18  is  a  diagram  of  the  acoustic  echosounder  layout. 

A  comparison  test  was  run  with  the  probe  system  and  the 
echosounder  between  13  00  and  203  0  hours  local  time  on  3 
September  1988  on  the  roof  of  Spanagel  Hall  on  the  grounds  of 
the  Naval  Postgraduate  School  at  Monterey  CA.  The  acoustic 
echosounder  was  placed  on  the  sixth  level,  northwest  corner 
of  the  roof  while  the  probe  was  attached  to  a  rigid  pole  and 
extended  approximately  1.5  meters  off  the  seventh  level  of  the 
northwest  corner  of  the  roof.  In  this  position  the  probe  was 
approximately  9.2  meters  above  the  echosounder  array,  thereby 
being  outside  of  the  echosounder  blind  zone.  Due  to  the 
building  itself  and  heating  exhaust  vents  on  the  eastern  side 
of  the  building  it  was  necessary  to  monitor  the  wind  direction 
to  insure  the  prevailing  winds  were  not  passing  over  the 
building  and  picking  up  heat  from  the  exhausts,  which  would 
have  greatly  affected  the  data.  Therefore  wind  speed  and 
direction  as  well  as  temperature  and  humidity  measurements 
were  taken  every  15  minutes  to  update  temperature  and  humidity 
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information  for  the  two  systems  as  well  as  determining  if  the 
prevailing  winds  were  flowing  across  the  building  before 
passing  over  the  instruments,  corrupting  the  data.  During  the 
entire  experiment  the  wind  shifted  several  times  but  it  was 
always  from  Northwest  to  Southwest.  At  all  times  the  air  flow 
passed  over  the  instruments  before  passing  over  the  building. 
All  the  data  was  valid. 

Figures  19  through  21  represent  a  portion  of  the  data 
collected  during  this  experiment.  In  each  of  the  figures  the 
upper  graph  is  the  echosounder  profile  of  the  atmosphere,  the 
central  graph  is  a  15  minute  time  averaged  CT2  profile  of  the 
atmosphere  based  on  the  data  collected  by  the  acoustic 
echosounder  and  the  bottom  graph  is  the  CT2  versus  time  plot 
measured  by  the  probe  system  at  an  altitude  of  9.2  meters 
above  the  echosounder.  The  dark  lines  below  about  6  meters  in 
the  upper  plot  and  the  discontinuities  below  6  meters  in  the 
center  plot  are  due  to  the  blind  zone  of  the  echosounder.  To 
compare  the  plots,  the  average  of  the  bottom  plot  was  compared 
with  the  value  at  an  altitude  of  9.2  meters  in  the  center 
plot.  Figure  19  was  taken  early  in  the  afternoon  and  shows 
strong  convective  pluming,  which  causes  a  higher  temperature 
structure  parameter.  Figure  2  0  which  was  taken  closer  to  the 
neutral  event  shows  a  marked  decrease  in  the  turbulence  and 
a  corresponding  decrease  in  the  values  for  CT2.   Figure  21 
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Figure  19 .  Echosounder  Trace  and  CT  Measurement  and 
Probe  Measurement  During  Strong  Turbulence 
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taken  during  the  neutral  event,  shows  virtually  no  turbulence 
and  a  much  lower  value  for  CT2.  Appendix  D  contains  additional 
measurements  taken  during  this  experiment  to  show  the 
corresponding  increases  and  decreases  in  CT2  for  both  the 
echosounder  and  the  probe  system.  The  purpose  of  these 
measurements  was  not  to  actually  sample  the  atmospheric 
processes  at  this  location  but  to  make  a  quantitative 
comparison  test  between  the  acoustic  echosounder  and  the 
probe  system.  For  a  complete  description  of  the  atmospheric 
turbulence  measurements  and  processes  for  this  location  see 
Weingartner  [Ref.  6]. 


B.   SCALE  SIZE  ERRORS 

The  temperature  variations  in  a  turbulent  atmosphere  range 
in  size  from  millimeters  to  hundreds  of  meters.  Optical 
aberrations  are  primarily  caused  by  variations  the  size  of  a 
Fresnel  zone  (AD)1/2  therefore  with  laser  frequencies  and  path 
lengths  of  several  kilometers  the  important  scale  sizes  are 
on  the  order  of  several  centimeters. [Ref .  13]  With  a  frequency 
response  of  150  Hz  and  an  average  wind  speed  of  2-5  m/sec,  the 
system  is  limited  to  scale  sizes  greater  than  3  cm,  which  will 
introduce  a  small  amount  of  "inner  scale"  error,  since  minimum 
scale  sizes  are  on  the  order  of  millimeters,  but  if  used  for 
measurements  in  conjunction  with  laser  propagation  through  the 
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atmosphere  the  error  will  be  negligible.  Additionally,  the 
acoustic  echosounder  utilizes  the  smaller  scale  sizes  of 
approximately  3  cm,  thus  there  will  be  negligible  error 
introduced  by  this  in  a  comparison  test. 

To  find  the  outer  scale  length  errors  we  can  express  the 
structure  function  of  the  probe  system  with  limiting  scale 
lengths  by, 


D(a,b)  =  4 


b 


(wrj)2    sin2(kr/2)k  5/3 


dk 


1+(UT1 ) 


l  +  (  UT  ^   ) 


(22) 


where  a  is  the  limiting  lower  frequency  (outer  scale) 

b  is  the  limiting  upper  frequency  (inner  scale) 

k  =  2n/A  where  A  is  the  actual  scale  length 

(jo   =   kV  where  V  is  the  wind  velocity 

r1  =  RC  time  constant  of  the  high  pass  filter 
(the  upper  frequency  cutoff) 

r2  =  frequency  response  of  the  probes 
(the  lower  frequency  cutoff) 

r  =  probe  separation  distance 

and  comparing  it  with  the  structure  function  over  all  scale 

lengths, 


D(0,co)  =  4 


.  2 

sin 

'  kr  ' 

0 

.      2  J 

-5/3 


dk  , 


(23) 


the  outer  scale  length  limiting  error  can  be  determined. 
Figure  22  is  a  graphical  representation  of  this  comparison, 
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showing  the  error  over  a  range  of  limiting  outer  scale 
lengths.  The  8%  limiting  error  on  the  low  end  of  Figure  22  is 
due  to  the  high  frequency  cutoff  of  the  circuit  and  as  the 
scale  size  decreases  the  larger  low  frequency  errors  of  the 
circuit  begin  to  dominate,  increasing  the  error.  The  limiting 
scale  length  for  the  experiment  can  be  determined  as  the 
height  above  ground,  which  was  approximately  30  meters.  From 
Figure  22  it  is  clear  the  error  introduced  due  to  finite  outer 
scale  lengths  is  approximately  12%.  The  design  of  the  acoustic 
echosounder  is  resistent  to  outer  scale  errors  therefore  the 
finite  inner  and  outer  scale  error  introduced  would  cause  the 
probe  system  to  record  measurements  approximately  12%  lower 
than  the  echosounder. 
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The  acoustic  echosounder  is  susceptible  to  inner  scale 
errors  which  will  cause  it  to  read  higher  than  the  probe 
system.  The  inner  scale  length  is  inversely  proportional  to 
the  wind  speed  and  can  be  expressed  as, 


1  =  7.4(„3/01/4.  (24) 


Ochs  and  Hill  [Ref.  29]  made  extensive  measurements  of  the 
inner  scale  length,  based  on  their  results  and  the  mean  wind 
speed  of  6  m/sec  during  the  measurements,  the  approximate 
inner  scale  length  was  3  mm.  At  the  edge  of  the  inner  scale 
of  turbulence  there  is  a  bump  in  the  temperature  spectrum  due 
to  diffusion  as  it  enters  the  viscous-convective  range.  Figure 
23  illustrates  this  bump  showing  the  spatial  power  spectrum 
$T  of  temperature  fluctuations  versus  the  scaled  wave  number 
Kri ,  which  is  the  wave  number  normalized  by  the  inner  scale 
length.  [Ref.  30]  Here  k  is  equal  to  27r/Scale  Length  and  r?  is 
equal  to  1^7.14  (for  air).  The  limiting  inner  scale  size  of 
the  echosounder  is  3.4  cm  (A/2  where  A=(340m/sec) * (5kHz) ) 
therefore  with  an  inner  scale  size  of  3mm  the  scaled  wave 
number  is  approximately  9xl0'2.  Figure  23  shows  the  acoustic 
echosounder  will  read  approximately  5%  higher  than  the 
Kolomogorov  spectrum  and  therefore  5%  higher  than  the  probe 
system. 
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Figure  23.  Spatial  power  spectrum  *T  of  temperature 
fluctuations  versus  scaled  wave  number  ktj  .    Solid  curve 
is  actual  model;  the  dashed  curve  is  Tatarski's  model. 
[Ref.  30] 


C.   ANALYSIS  OF  DATA 

Figure  24  is  a  comparison  of  15  minute  time  averaged  data 
collected  from  the  echosounder  and  the  probe  system.  This  data 
was  taken  before  the  noise  measurements  and  discovery  of  the 
5  millivolt  DC  offset  error  in  the  A  to  D  converter  RMS  module 
combination.  Figure  25  shows  the  corrected  data  comparing  the 
two  systems.  The  data  clearly  shows  the  correlation  of  the  two 
systems  even  with  the  volatile  trends  of  the  turbulent 
fluctuations.  It  also  shows  a  decrease  in  the  temperature 
structure  parameter  leading  up  to  and  during  the  neutral 
event,  which  corresponds  with  the  actual  physical  processes 
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Correction  for  DC  Offset. 
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Figure  25.  Comparison  of  Data  From  Acoustic 

Echosounder  and  Temperature  Probe  After 
Correction  for  DC  Offset 
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TABLE  3 


MEASUREMENTS 


(CORRECTED  FOR  DC  OFFSET) 


TIME (LOCAL) 

ECHOSOUNDER 

PROBE 

%  DIFFERENCE 

1400 

.0354184 

.0368390 

3 

1415 

.0316086 

.0391012 

23 

1430 

.0358514 

.0335893 

6 

1445 

.0500970 

.0492531 

2 

1500 

.0322137 

.0321643 

0 

1515 

.0490100 

.0432243 

12 

1530 

.0474621 

.0384319 

20 

1645 

.0311984 

.0271304 

13 

1700 

.0263566 

.0350252 

32 

1715 

.0361568 

.0292298 

20 

1730 

.0260007 

.0232678 

10 

1745 

.0217524 

.0274032 

25 

1800 

.0233976 

.0200511 

14 

1815 

.0141386 

.0165680 

17 

1830 

.0138539 

.0165680 

18 

1845 

.0132466 

.0184468 

39 

1900 

.0129296 

.0145711 

15 

2000 

.0054777 

.0077345 

40 

2015 

.0058262 

.0071837 

23 

AVERAGE  %  DIFFERENCE 

17 

going  on  at  the  time.  With  the  approach  of  sunset,  at  193  3 
local  time,  the  sun  heated  the  earth's  surface  to  a  lesser 
degree  thereby  reducing  the  temperature  difference  between  the 
earth's  surface  and  the  air,  which  in  turn  reduced  the 
temperature  structure  parameter  and  the  turbulence. 
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Table  3  contains  the  values  of  CT2  measured  by  each  of  the 
devices  and  corrected  for  the  offset,  it  indicates  an  average 
difference  of  17%,  with  the  probe  system  reading  lower.  The 
scale  length  errors  indicate  the  probe  should  read 
approximately  12%  lower  due  to  outer  scale  length  errors  and 
5%  lower  due  to  inner  scale  bump  errors.  Therefore  there  is 
no  significant  difference  between  the  probe  system  readings 
and  the  echosounder  calculations. 
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V.  CONCLUSIONS 

Independent  verification  of  CT2  values  measured  by  the 
acoustic  echosounder  is  important  [Ref.  28]  and  the 
differential  temperature  structure  parameter  probe  has 
provided  a  valuable  comparison  indicating  the  absolute  CT2 
values  of  both  the  echosounder  and  the  probe  are  valid.  Taking 
into  account  all  known  errors  there  is  no  significant 
difference  between  readings  of  the  echosounder  and  the  probe 
system,  which  is  an  extremely  good  indication  that  both 
systems  are  providing  valid  measurements.  Additionally  this 
thesis  demonstrated  that  solar  heating  of  the  probes  in  the 
atmosphere  does  not  appear  to  play  as  significant  a  role  as 
first  thought.  The  only  major  effect  solar  heating  has  on  the 
differential  system  is  when  one  probe  is  directly  illuminated 
by  solar  radiation  while  the  other  is  shaded. 

This  probe  system  has  many  applications  including  being 
placed  on  towers  to  calibrate  other  turbulence  measuring 
devices  as  well  as  being  attached  to  a  rawinsonde  system  and 
launched  to  measure  the  vertical  profile  of  CT2.  When  used  in 
this  mode  it  can  measure  values  of  CT2  up  to  30  km  altitude 
accurately,  however  if  it  rotates  as  it  ascends  through  the 
atmosphere  the   sun/shade  effect   of   solar  heating  will 
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adversely  affect  the  system.  Further  developments  to  the 
package,  such  as  addition  of  wind  vanes  on  the  probe  assembly 
which  will  not  affect  the  turbulent  flow  but  will  dampen  the 
rotation,  or  a  small  motor  with  a  flywheel  to  act  as  a 
gyroscopic  stabilizer  to  prevent  the  probe  from  rotating,  will 
eliminate  the  errors  induced  by  this  effect. 

Other  improvements  to  the  system  include  methods  to 
automatically  update  the  temperature  into  the  program  from 
the  balloon  systems  onboard  temperature  sensor.  Another 
improvement  would  be  to  increase  the  data  transmission  rate 
of  the  rawinsonde  telemetry  system  to  get  a  higher  resolution 
profile  of  the  thin  stratified  layers  of  the  turbulent 
atmosphere  or  even  possibly  having  the  system  transmit  an  AC 
signal  from  which  a  great  deal  more  information  can  be 
extracted  such  as  the  power  spectral  density  of  the 
turbulence. 
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APPENDIX  A 

PROBE  SYSTEM  SOFTWARE 

The  program  that  runs  the  system  is  called  "CTSQR".  It 
controls  the  probe  system,  collects  the  data,  reduces  it,  and 
then  displays  and  stores  it  for  further  analysis.  It  is  based 
on  the  same  program  that  controls  the  acoustic  echosounder. 
It  can  be  broken  down  into  several  sections.  The  first  section 
sets  up  the  system,  initializes  all  arrays,  creates  a  data 
file  which  can  store  up  to  eight  hours  of  data  and  sets  up  the 
function  keys  which  are  used  to  update  the  temperature  used 
in  calculating  the  Seebeck  Coefficient,  prints  out  raw  voltage 
data  or  ends  the  program  storing  what  has  been  collected.  The 
next  section  initializes  the  Infotek  AD200  analog-to-digital 
converter,  which  collects  the  data.  Now  that  the  system  is 
ready  to  collect  data  it  calculates  the  Seebeck  Coefficient 
based  on  the  information  input  at  start-up  or  updated  through 
the  function  key.  Next  it  collects  data  every  second  and 
averages  it  over  ten  seconds,  reduces  it  to  CT2  and  plots  it 
every  ten  seconds.  Every  15  minutes  it  prints  out  the  plot  and 
then  resumes  the  data  collection.  The  program  "CTREADER"  can 
take  the  data  file  generated  by  "CTSQR"  and  read  it  and 
calculate  15  minute  time  averaged  values  of  CT2. 
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19  I          RE-STORE  'Ctsqrz  .700,l,e• 
20  i  ctsor  is  sep  isee-  mro 

39  I  This  preoren  collects  one  date  channel  fron  a  HP  3421ft  or  ftO  converter, 

4?  I  end  stores  eight  hours  of  the  binary  data  on  a  disc  file 

5?  • 

G?  OPTION  BftSE  1 

7?  I  Initialize  the  arrays 

PP  DIM  Pes»M6J,Disc_address*t.20J,Fllel*r30) 

9?  INTEGER  I.J.J4,K,Kstert ,Kend.N,Nrec,Hr.O2(2BB0,4),Plolnun,Prlnt_key 

1??  I   02  •  The  reduced  data  output  array  (28B0.4) 

11?  i      -  (Pay.Hr,  VOLTS_ftUS,  CT  ) 

12?  INPUT  "ENTER  ftIR  TEMP  7DE6REES  C)",T 

13?  ' 

14?  I   Set  constants 

15?  Dl8c_addresst-'- ,700,e,6"  I  HPIB  address  of  disc 

IE?  Galn-500??                I  ftnpllfler  gain 

17?  Ma-»rec-2ee?               I  t  records  In  output  file 

ie?  Nplot-9??                I  *  points  plotted 

19?  Plotnun-0 

2"?  Print_key-1               I  print  raw  data  If  >0 

21?  R-l.?                    I  Probe  separation  (n) 

22?  Scale-1????          I  Disc  storscje  scale  factor 

2"?  I  The  Equation  For  The  Seebeck  Coefficient 

21?  See-3.B7?7E-2»e.534BE-5«T-3.3l35E-7»T"2 

25?  Beck--2.77432E-9'T-3-1.2S3E-ll«T*4 

2E?  Seebeck-(See+Beck >«1 .E-3 

27?  PRINT  Seebeck 

2e?  I 

29?  R_one_lhlrd-R'( 1./3.  ) 

3??  I 

31?  "SET  TIMEDftTE 

32?  INPUT  "00  YOU  WANT  TO  RESET  THE  CLOCK  (Y  OR  N>7",0« 

33?  IF  0*-'Y-  THEN 

319  INPUT  "ENTER  "DO  MMM  YYYY"  (Local  Tine  )', Date* 

25?  INPUT  'ENTER  "HR-MIN  SC"  (Local  Tine  )',  Tine* 

35?  SET  TIMEOATE  0ATE< Date*  >+TIME( Tine* ) 

27?  PRINT  DATE*(T1MEDATE>,TIME*(TIMEDATE> 

30?  Tstart-TIMEPATE 

3D?  T?-Tstarl  MOO  8B409 

40?  FND  IF 

11?  INPUT  "INPUT  SITE  NAME", Site* 

12?  INPUT  -ENTER  THE  ft-D  CONVERTER  0FFSET(-.?95  FOR  HP  2l7)",Zero 

41?  ItlPUT  "ENTER  THE  LOWEST  DECftDE  FOR  THE  PLOT  ( NORMAL  USE  -Sl'.Ynln 

11?  i 

4^9  Crp«lr_f I le     I 

4RP  '  Set  up  the  dMa  reduction  output  file 

17?  INPUT  "ENTER  REOUCEO  DftTft  OUTPUT  FILE  NAME'.FlleU 

IP?  Fllel*-Fllel*SDiec_address* 

4  "?  i 

5?0  INPUT  "1ST  ENTRY  IN  REDUCED  OUTPUT  FILE7  (YES  OR  N0)-,0» 

51?  IF  0*-'N0-  THEN  GOTO  Oldflle 

52?  Neufile'  CREftTE  BOAT  Fl lel», 1 , 23949     I  2  BYTES  *  FILE  SIZE  B  HOURS  OF  DATA 

530  ASSIGN  BFllel  TO  Fllel* 

519  Nrec-9  II  OF  ENTRIES  IN  THE  OUTPUT  FILE 

SB?  GOTO  Setup 

SB?  Oldflle^  ASSIGN  PFUel  TO  FllelS 

570  ENTER  PFllel ;02(«  ) 

5P0  ASSIGN  "Fllel  TO  Fllel*        IGO  TO  START  OF  FILE 

590  Nrec-D2(l,l)  1  THE  OLD  *  OF  ENTRIES  IN  THE  OUTPUT  FILE 

BP9  Setup   i   Set  up  the  dote  reduction  and  plot  fornet 

G10  OUTPUT  KBD; 'SCRATCH  KEYE";          I  Clear  keys 

G20  CONTROL  2,2; I                       I  Select  user  nenu  1 

G30  ON  KEY  1  LABEL  "PRINT     RAW  60T0  Prlnt_raw 

G4?  ON  KEY  B  LAPEL  -Quit"  GOTO  Quit 

GT.?  ON  KEY  2  LAPEL  "UPOftTE   TEMP*  GOTO  Upnate_lenp 

G5?  ' 

G7?  Npolnt-TIMEDATE  MOO  8G400  MOD  3600  DIU  15 

E00  CALL  Plot5etup(Nplot,Slte*,Ynln) 

G90  CALL  Inlt_ad20?                     1  Initialize  A-D 

700  •  OUTPUT  KBD;  "L-;                     I  Turn  on  graphics 

710  I 

720  '  Begin  the  naln  data  acquisition  loop 

730  • 

740  WHILE  NreOMaxrec 

750  Start  lo    I 
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76?  Voltso-0 

770  Store_data-0 

790  FOR  1-1  TO  10 

790  Sync;     I  Synchronize  data  collection  with  the  syeten  clock 

800  Tl-INTUIMEDATE   MOO   96«00> 

810  IF    TKT0   THEN   T0-T0-B6400 

820  IF  (T1-T0X1  THEN  GOTO  Sync 

830  T0-T1 

840  Reed_ed   I 

850  CALL  Adln(Voltage,Zero>           I  Read  Infotek  A-D 

860  IF  Prlnt_key>0  THEN 

870  PRINT  USING  ' 1B0.D0DD' ; Vol teoe 

880  GOTO  Read_ad 

890  END  IF 

900  Vol  t  sq«Vol  tsq+Vol  tage'Vol  tega    I  Average  voltage'2 

910  Npolnt-INTITl  HOD  3600  MOO  Nplot  ) 

920  IF  Mpolnt<Npotnt_old  THEN  Store.data-l 

930  IF  Npoint_old>0  AND  Store_dala-0  THEN 

9*0  I  Plot  the  data 

950  CtsQr-<Voltege/(6aln«Seebeck«R_one_thlrd)>*2' 

960  Lgctsqr-LGT(Ctsqr) 

970  l  MOVE  Npolnt_old,Volts_old 

980  l  DRAW  Npolnt, Voltage         I  Plot  Voltage 

990  MOVE  Npo!nt_old,LgctsQr_old 

1000  DRAW  Npolnt, Lgcteqr 

1010  END  IF 

1020  Velt9_old-Voltage 

1030  Lgctsqr_old-Lgctsqr 

1040  Mpolnt_old-Npolnt 

1050  IF  Store_data-l  THEN  Npolnt_old-0 

1060  NEXT  I 

1070  Vclts_avg-SOR(Voltsq/10) 

10e0  Ct-Volts_avg/<6aln«Seebeck»R_cne_thlrd> 

1090  DISP  •  ' 

1100  Tl-TIMEPATE 

1110  Dey*-PATE*(T1  ) 

1120  Tlne*-TIME*f Tl > 

1130  • 

1  140  l  •••  CT  SQUARED  DATA  REDUCTION  SECTION  ••• 

1150  I 

1160  Yrt-Deytie, 11) 

1170  Day-((DATE(Day»)-DATE( '1  JAN  '&Yr«>)  DIV  BG400J+1 

lie0  l 

1190  Nrec-Nrec+1 

1200  '  CALCULATE  OECIMAL  HOURS 

1210  TS-Tlne* 

1220  Hours -VAL<T*(  1,2)  >♦<  VAL<  T»t  4,5)  HVAL  <  T«I  7,8  )  )/60  )/60 

1230  Hr-1000»Heurs  !  Note  that  the  HP  rounds 

1210  PRINT  'Record  f;Nrec:  "  Collected  "iDay*-;-  "iTlne* 

1 250  PRINT  ■  • 

1260  ALPHA  OFF 

1 270  I  SET  UP  OUTPUT  ARRAY 

I2e0  Nl-Nrec+1  IFOR  OPTION  BASE  1 

1290  D2( 1 , 1  )"Nrec 

1300  D2(Hl,l)-Day 

1310  D2UI1  ,2>-Hr 

1320  D2<N1.3)-Volts_avg«Scale 

1330  D2(H! ,4)-Ct»Sca!e        '  Intensity 

1340  PRINT  Nl ;D2(N1,1  )  ;D2( Nl , 2  ) : D2( Ml ,3  )  ;D2( Ml . 4  ) 

1350  Stcre_dBta  !  Write  output  every  NPLOT  seconds 

1360  IF  Stcre_data-1  THEN 

1370  Store_dele-0 

J3P0  PRINT 

1390  PIS"  "WRITING  REDUCED  OUTPUT' 

1400  OUTPUT  fFllel ;D2( O 

1410  ASSIGN  PFllel  TO  FllelS 

1420  Plclnun-Plo(nun*l 

1470  IF  Plctnun  MOD  2-1  THEN 

1440  PRINTER  IS  701 

1450  PRINT  • 
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I  Dump  screen  to  printer 


1460  PRINTER  IS  1 

1470  ENO  IF 

ltee  DUMP  GRAPHICS  1781 

1490  Npclnt-0 

1500  PISP  '  ' 

1510  CALL  Plotsetup<Nplot ,Slte«> 

1520  EMD  IF 

1530  End_uhlle   END  WHILE 

1540  I 

155?  Updete_tenp-   •  Updates  Seebeck  Coefficient  With  New  Air  Tenp 

1560  INPUT  "ENTER  NEW  AIR  TEMPERATURE? DEGREES  O'.T 

1570  See-3.B707E-2+8.534BE-S»T-3.313SE-7«T*2 

ISP?  Beck— 2.77432E-9'T-3-1.253E-U«T-4 

1590  Seebeck-(See+Beck  JM.E-3 

IBM  GOTO  Start_lo 

1610  Prlnt_rew  I  TOGGLE  THE  PRINT  FLA6 

1E20  Prlnt_key— Prlnt_key 

163*  GOTO  Star-t_lo 

1 E if  Quit    I  Urlte  reduced  data  output  file 

JES0  FOR  1-1  TO  Nl 

1660  PRINT  I:D2(  1.1  );D2(  I,2);D2<  1 ,3  )  ;D2(  I ,  *  )  ;D2(  I,5);02(  1,6) 

1G70  NEXT  I 

16e0  OUTPUT  PFllel;D2(» ) 

1690  PRINT  "DATA  FILE  HAS  BEEN  STORED  UNDER  NAME'.Fllel* 

wee  BEEP 

1710  BEEP 

1720  ASSIGN  BFllel  TO  t 

1730  STOP 

1740  END 

1750  SUB  Plctsetup<Nplot, SlteS, Ynln) 

17E0  Ynax-1 

1770  GIMIT 

]7e0  GRAPHICS  ON 

1790  LINE  TYPE  1 

ie00  VIEWPORT  15,120,15,60 

iei0  WINDOW  0, Nplot, Yfiln.Ynax 

ie20  AYES  B0,.5.0,Ynln,5,2 

1830  CLIP  OFF 

IP40  CSIZE  4, .6 

ieS0  LORG  G 

ieG0  •    Draw  Log  Y  Axis 

ie70  I    FOP  Pecade-Ynln  TO  Ynax 

leee  <       FOR  Unlts-1  TO  H8»<Decade<Ynex) 

ie90  i       Y-Oecade+LGT<Unlts) 

1900  i       MOVE  0.Y 

1910  '       DRAW  Nplot ,Y 

1920  i       NEXT  Units 

1930  i    NEXT  Decade 

1940  i   Label  horizontal  axis 

1950  T)=T1MEDATE  MOD  BG400 

19S0  Hrs-Tl  DIV  3E00 

1970  T2-T1  MOD  3B00 

1900  Mln-T2  DIV  G0 

1990  Qtrhr-Mln  DIV  IS 

2000  FOR  M-0  TO  Nplot  STEP  300 

2010  MOVE  M,Ynln-.06 

2020  Otrnln-Otrhr»15+(M/300)»5 

2030  IF  Otrmln-60  THEN 

2040  Qtrmln-0 

2050  Hrs-Hre+i 

20B0  END  IF 

2070  LAPFL  USING  "DD, A, ZZ" :Hrs :"•• ;Qtrnln 

20P0  NEXT  M 

2090  MOVE  Nplot/2,Ynln-.B 

2100  LABEL  "Tine  (Local  )" 

2110  '  Label  Ordinate 

212e  LORG  8 

2130  FOR  M-Ynln  TO  Ymex 

2140  LORG  9 

2150  CSIZE  4 

21E0  MOVE  -Nplot/23.3,M 

2170  LABEL  USIN6  't,Y.':m\Z' 

2ie0  CSIZE  2 

2190  LORG  1 

2200  MOVE  -Nplct/2e,M 

2210  LABEL  USING  '!,*■  ;M 

2220  NEXT  M 
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rVL'tf  NfcXI    M 

:t>ne  hove  Npiot/2,Ynm-.o 

2100  LAPEL  "Tine  (Local  )' 

2110  I  Label  Ordinate 

2120  LORG  B 

2131  FOR  M-Ynln  TO  Yna* 

2140  LORG  9 

2150  CSIZE  4 

2160  MOVE  -Nplot/23.3,M 

2170  LABEL  USING  *t,K*:*l«* 

2JP0  CSIZE  2 

2190  LORG  1 

2200  MO'.'E  -Np)ct/28,M 

2210  LABEL  USIN6  "t,K";M 

2220  NEXT  M 

2230  LD1R  PI/2 

2240  LORG  G 

2250  CSIZE  4 

22G0  MOVE  -Nplot/7,<Ynln+Y«ax  )/2 

2270  LABEL  "L06  OF  CT"2" 

22B0  ITltle  the  plot 

2290  LOIR  0 

2300  LORG  4 

2310  MOVE  Mplct/2,Ynax+l 

2320  LABEL  DATEK  TIMEOATE  ) :  "     ":Slte* 

2330  CLIP  ON 

2340  SL'BENO 

2350  SUB  AdlnfJcltage.Zero) 

23G0  !  2G  APR  19PG'  OLU 

2370  I  IMFOTEK  A-0  Input  routine   set  up  for  Internal  trigger 

2360  !  and  average  40  points  over  three  G0  Hi  cycles 

2390  ! 

2400  INTEGER  I .Npolnt s, Ad_data< 1 • 40 > 

2410  DIM  Select*[20) 

2420  Ad_sel_ccde-17 

2430  lntenslly-0 

2440  Mpoints-40 

2450  Count*-VAL»<Npolnts) 

2460  Scale-5.0/(Npolnts«2047.  ) 

2470  Stdev-0 

2<e0  Del ta_tlne*-" 1250000"       1  AO  Interval  between  sanples  In  nsec 

2490  Select*-"select  Isl  end" 

2500  GOSUB  Read_ad 

2510  FOR  1-1  TO  Npolnts 

2520  Vol tege-Vollage+Ad_data( I > 

2530  NEXT  I 

2540  Voltage-! Volt age- Zero) 'Scale 

2550  GOTO  Subend 

2560  Read_ad      I  Read  the  Infotek  A-D 

2570  i  Initialize  the  A-D 

25e0  OUTPUT  Ad_sel_code: "RESET", "Internal", "count", Count* 

2590  OUTPUT  Ad_sel_code;"tine".Delta_tlne»,"deleyon", Select! 

2600  OUTPUT  Ad_sel_code: "STATUS" 

2G10  ENTER  Ad_sel_code ;Resp» 

2G20  IF  Resp»-" "  THEN 

2630  ENTER  Ad_sel_code  USING  "»,U" ; Ad_data< •  ) 

2640  OUTPUT  Ad_sel_code; "STATUS" 

2650  ENTER  Ad_sel_code : Respi 

2GG0  IF  Resp»<>" "  THEN 

2G70  PRINT  "ERROR  DURING  SAMPLIN6  -  ";Resp« 

26P0  END  IF 

2690  ELSE 

2700  PRINT  "ERROR  DURING  A-D  INITIALIZATION  -  ";Resp» 

2710  END  IF 

2720  RETURN 
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2730 

Subend    1 

2740 

SUB END 

2750 

SUB  Inlt_ad200 

27S0 

UNITIALIZE  AO_200 

2770 

Code-17 

2780 

Dunny-REA0I0( Code , 3 ) 

2790 

URITEIO  Code, 010 

2800 

CONTROL  Code,0;l 

2810 

SUBENO 
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APPENDIX  B 

SOLAR  HEATING  PROGRAM 
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RE-STORE  "ALTPL0T1  ,738,1,9" 

i  "»  nrj     i  goo   f*RQ 

THIS  PRCSRAM  PLOTS  THE  CHAN6E  IN  THE  DELTA  T  '.  TEMPERATURE  DIFFERENCE 
BETWEEN  THE  THERMOCOUPLE  AND  THE  AIR  TEMP  CLE  TO  SO'  A.9  H£ATINS>> 
WITH  INCREASE  IN  ALTITUDE  2ASE0  ON  SCCC'S  CALCULATIONS  IN  AFSL  P"2 
DATED  22  FEB  1284  AND  CAMPBELL'S  WORK  FROM  OCT  I2S9 


i  THE  FCLLCUIMS  IS  THE  LIST  OF  VARIABLES 
i 

J "4 . IS 

S-.l* 

Ep»i-.25 

E3S3-.2S 

Fd-PI 

Fr-2. 

Ep»l»-.S 

Esai-.S 

Rs-.322 

Ra-.89S 

Rg-.315 

Albedo-. 25 

Sigm>s-5.S7E-12 


CM-. 8883254 

0cw.ee254 


Beta-1.4S2E-S 

Suth-118.4 

Pr-. 714 

Re   « 

Nud  - 

T 

Rex  - 

Mu   - 

Nu   - 

K 

H 


Alt  - 
V  — 
Kr  - 
M 

Oelt- 
Deltl- 


mechanical  equivalent  at   heat  (W  Cal*-1  Sm ) 

Solar  Constant  (W  cm*-2 > 

Visible  Wavelength  ABSORPTIVITY  of  th*  THEPMOCCL'PLE 

Short  Wave  Emisaivity  for  CAMPBELL'S  Calculations 

Forr<»  Factor  for  Direct  Radiation 

Fern  Factor  for  Reflected  Radiation 

Long  Wave  EMISSIVITY  of  the  THERMOCOUPLE 

Long  Uave-  Emissivity  for  CAMPBELL'S  Calculations 

Short  Wave-  incoming  Radiation  (Cal  crT-2  see*-l> 

Long  Uave  Atmospheric  Radiation  (Cal  cm"-2  sec*-l> 

Long  Wave  Radiation  from  Ground  (Cal  em"-2  sec*-l> 

Reflection  ef  the  Earth 

Stefan-Bo Itsmann  Constant  (W  e»*-2  K*-4> 


Slgeal-1.2SE39E-12  \    Stefan-Bolt  :wnn  Constant  (Cal  cm*-2  K*-4) 


Prebe  Diameter  in  meters 
Probe  Diameter  in  centimeters 


Ccnvert-4. 12992Ei2  I  Conversion  Factor  for  Thermal  Ccndus 


Constant  for  Determining  Mu  <kg"-l  n"-l  K*-l/2> 

Sutherland's   Constant    for  Mu   <K> 

Olmensionless  Prandtl  Number  fcr  A;r 

Oimensioniess  Reynolds  Number  for  air 

Oimensionless  Mussel t  Number  for  Air 

Air   Temperature   (Kelvins) 

Air  Density  as  Function  of  Altitude  (kg.'m"2> 

Dynamic  Viscosity  of  Air  (N-see/m-2 > 

Kinematic  Viscosity  of  Air  (m'2.'sec> 

Thermal  Conductivity  (Cal  sec"-!  cm"-2  (C'cmr 

Convectl.e  Heat    Transfer  Coefficient    tfiALTM 

(Cal    em*-2   sec*-l    K"-l > 

Percentage  of  Solar  Radiation  Reaching  6lven  Altitude 

Heat  Flux  from  Earth  at  Siren  Altitude  (W  cm*-2> 

The  Slven  Altitude  <Km> 

Wind  Speed  (m.'s*c.v 

Empirical  constant  based  en   Reynolds  Number  from  Krleth 

Empirical  constant  based  on  Reynolds  Number  from  Krleth 

Temperature  Difference  between  TC  and  Air  for  SOCD 


■1 


Temperature  Difference  between    TC   ar.c 
Both   m  Degrees  C 


'SELL 


INPUT  "WHAT  WIND  SPEED  00  YOU  WANT'.V 
l 

•  PLOT  SET'-'P 

INPUT  "INPUT  1  FCR  PLOTTER  OR  2  FOR  CRT",C 

IF  0-1  THEN  PLOTTER  IS  787,'HPSL" 
t  BE  SURE  PLOTTER  OIP  SWITCHES  ARE  PROPERLY  SET  is 

IF  0-2  THEN  PLOTTER  IS  CRT ," INTERNAL" 
i  SINIT 
'  GRAPHICS  CM 

LINE  TYPE  1 

VIE'JPCRT  15,128,15.98 


SWITCH  1,2.2  IN  POSIT  1 
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553  UINCOU  2,22,8, 1 

573  AXES  5,. 35, 2, 2, 5, 2 

cog  CLIP  OFF 

com  f-c  T  7C  4  _  c 

SC2  LCRS  S 

512  !  LABEL  HORIZONTAL  AXIS 
523  FOR  M-0  TO  23  STEP  5 
523  MOVE  M.-.84 

5*0  LABEL  M 

553  NEXT  M 

553  MOVE  15  - . 1 

573  LASEL  "ALTITUDE  <Km!" 

582  '  LASEL  VERTICAL  AXIS 

593  LOPS  9 

732  FOR  M~9  TO  1  STEP  .  ! 

713  MOVE  —  .5  M 

722  LASEL  M 

723  NEXT  M 
7*0  LOIR  PI/2 
753  L0R6-  S 

7C2  MOVE  -4.5  .5 

772  LABEL  'DELTA  T  !C!" 

793  '  TITLE  PLOT 

792  LOIR  2                                             " 

939  LCRS  4 

910  MOVE  15  1 . 1 

923  LASEL  "SOLAR  HEAT  INS-  OF  THERMOCOUPLE" 

923  MOVE  15,1 

940  CSI2E  2 

953  LASEL  "FOR  A  WIND  SPEED  OF"  ;'„';  "m/see" 

953  CSI2E  4 

973  !• 

S90  \      Calculations 

993  ! 

903  Aiteid~9 

513  Oelteld~0 
923  Deitlcld-9 

S23  FOR  Alt-0  TO  23 

940  ►  I  an  assuming  the-  relat  lonship  far  F  and  Se  ars   linear- urt  Altitude 

950  «  The  Values  were  taken  from  BROWN  and   6C0D 

9S0  F~.S+<Alt*. 312902! 

973  Se-.BS6-<Alt*. 00141935! 

990  (  The  relationships  for  T.Rcw.Mu.Nu,  and  K  are  taken  from  the 

990  !  HANDBOOK  of  GEOPHYSICS  and  the  SPACE  ENVIRONMENT  chap  14 

1803  !•  and  based  on  a  U.S.  Standard  Atmosphere 

1013  •  This  is  taken  frc-m  the  Standard  Atmosphere  Temperature  Profile 

1023  IF  Alt<-10.  THEN 

1320  T-289.15-7.81S»Alt 

1843  ENO  IF 

1050  IF  Alt  MB.  AND  Alt<2»  THEM 

1059  T-219. 

1070  END  IF 

1890  IF  Alt>-23.  THEN 

1B98  T-219.+1.2MAlt-20. > 

1133  ENC  IF 

1113  Rcu~l  .22S2-.3299»Alt  '.    Density  Charge  is  apprc:<  inately  Linear  art 

1120  \                                                  Altitude  up  to  about  1C0  Km 

1120  Mu-(Seta»T-l.S!.'(T+Suth> 

1140  Nu-»Mu/Rsw)». 1224255     1^  Conversion  faster  to  get  correct  viscosity 

1150    Powei 12. /T 

1153  K1-2.S5219S-2»T-1  .5 

1172  K2-T+24S.4«10-Power 

1190  K-(K1/K2  !. 'Convert 

1193  I  The  determination  of  H  »as  based  on  KRAMERS  (1945!  which  was  shewn 

1220  I  to  mere  closely  approximate  experimental  results  by  the  plot  KPL0T2 

1213  Re-Rew»V*0,-i,'Mu 
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1222 
12*3 
1252 
'■"id 

23 

■ma 

i2:e 

1223 

i  mm 
1243 

I  7Cfl 

■Jen 

.272 
L253 
1292 
[422 
1412 
1 422 
[422 
.443 
[452 
1452 
i  47m 

1422 
[492 
1  e*>8 
12 
[522 
1522 
[E42 
1552 
L5S2 
1573 

[c92 
t  c2°* 
I c  '  2 
[S23 
1522 
1942 
[553 
[££3 
[£"2 
1  com 

[£92 
1  ""22 


'  7C2 
[752 
72 
[79C 
i  -'93 
1  922 
ISifl 
[9*'2 
[222 


ic  9«^cg  THEM 

Nud~. 21 »Pr"  .31»Rs*  .225 
EL£E 

Nud-.S»Pr*.21*Rs".E 
END  IF 
H-Nud»K/Dcn 
Calculaticna  baas-d  en  BRCUN  and  500D  (AFSL  Pub  1924! 

PRINT  \ALT";Alt  .•RCy";Rsi»,"M'J";Mu,*NU"  iNu 

PRINT  '    K";K,"Re-";Rs, "H";H, "T";T 

PRINT  "KRIETH'S  H-" ; CI ,  "FIRSTPART~" ; Al , "2ND-" ; Bl 

B-F*S»Ep»t»< !  1  ./FdHlAlbedc/Fr  ! ! 

C-£p3lw»«  (Es.'Fr  !-;£i5na»TA4!  > 

Ds-1 t-A*<B+C !  * 

LINE  TYPE  9 

MnMC  Altaic  Caltc-id 

DRAW   Alt .Dsii 

Calculations  baaisd   en  CAMPBELL    !19S9! 


I    Krl*th'a~«npirieal    farrt  of    1/h 


IF  Rs-<4  THEN 

KP-.S91 

N-.222 
END  IF 
IF  Rs<43  AND  Ra->-4  THEN 

u  r  «  o  7 ! 

KJ,  TOC 

END  IF 

IF  Ra->-42  THEN 
Kr-.SIS 
N-.465 

END  IF 
">.1-E333*R3*(  l  +  <  (PI'Albsdc  >.'2  > ) 

Bl^PI»Ep3l»<<(Ra+RB>/2!-(Sla=al»TA4 ) > 

Cl-iDcn.MKr'K  !  >-<  Nu/(  '„,»D«>  >"N 

Col  tlH  Al+Bl  )*C1 

LINE  TYPE  4 

MOVE  Aliold,D»ltlold 

DRAW  Alt ,D»itl 

Altold-ftli 

Ds-ltloid-Os-ltl 

Dsl tcld-Dsii 

PRINTER  IS  721 

PRINT  'ALT-'  i  Alt ,  "DELTC  AFSL  >-"  i  Ds-1 1 .  "0ELT(  CAMPBELL  >»*  i  Dsl  1 1 
NEXT  Alt 
rene  1 
1  Ti.jc  yvpc  g 
MOVE  1..94 
"RAW  2     .  94 
CRAW  2, .94 
DRAM  4  ,  . 94 
MOVE  7   94 
LINE  TYPE  1 
LABEL  "AFGL  FORK" 
LINE  TYPE  4 
MO'.'E  IS  .24 
DRAW  21  ,  .94 
LINE  TYPE  I 
MOVE  25,. 94 
LABEL  "CAMPBELL  FORM" 
DL'MP  GRAPHICS  *731 
EMC 
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APPENDIX   C 

CONDUCTANCE    PROGRAM 

i  n,        i  DC.CTnqc     "UOi  pt-j  .     nnm     i     nf 

_  y  !  t  *_.     w  .  wH_.         ",_,_,_.       t,wv.A,- 

""*  '    THIS    PROSRAM   COMPARES    THE   RESULTS   OF    THE   COMPUTATION   OF   THC   coj»j"cr»T t i ic 

tn  i     ppkjni  ip  Tiwrr    ov    c    .   raMOOci  i    'c    mctuoP    Gkin    qv    tuc    accj      o&pco'c    MCTunn 

^  nt  i     akin    oi  nrc     tucm    _iC_iTk|CT    cvproiMCMTii      iiai  iicc    cdpm    rAMppci  i 

cm  I      <  ">     nrT     1  ooo 

WW  A   A.  WW    I  A     WWW 

cm  i  iiai  iicc  aoc  cod  a  crn  aTMncouCDC  i  ._   1  ATM  at  ->fl»  ncc  r 

_,  W  ViiLwLJ        lini.        I      wl  I        li        W    I    — '        til    iiuji      lltllL         4.W*  1         lilil        III  kll        L.  _.  W        W 

"7171  I 

23  Row^l.22     i!<5/n"3   AIR  OENSITY 

S3  Mu-1.7.E-E   '  (N-sae >/n"2  CYNAMIC  VISCOSITY  OF  AIR 

133  Nu-1 .7SSE-S  !  n*2/3  KINEMATIC  VISCOSITY  OF  AIR 

113  K-S.175E-5   '        THERMAL  CONDUCTIVITY  OF  AIR 

123  D--2.54E-3    '  n  DIAMETER  OF  THE  PR08E 

123  Pr^. 714      !  OIMENSIONLESS  PRANDTL  NUMBER  OF  AIR  AT  23  DES  C 

143  '  KK  AND  N  ARE  EMPIRACAL  NUMBERS  BASED  ON  THE  REYNOLDS  NUMBER 

tern  i   cc    pccpdtocp  hi  i/otctu  iocc  —  <i  i 

A    WW  I   IW         W_.WWI«AW_-W  All         IllltC    IN  A    WW  W         W   W.         +U 

ICJ  I   II  TC  tuc  uTMr  coccn  T  kj  «./- 

»    W  \J  If  AW  I    I    IL         «    1'lt  W  I      *_C  W  ill  '"      J 

1-701  i     d-s.    tc     tuc    nTMCMCTnuji  pec    ocvmpi  nc    miimocd 

A     IV  I  I  _T        AW  i   IIL.        WAII_,l"WAW1«_.CwW        ii_.itlw_.ijw       IIUHUUII 

123  '  Nud  IS  THE  OIMENSIONLESS  NUSSELT  NUMBER  AS  DESCRIBED  IN  KRAMERS  194S 

193  Ccrrac.^.7E2413822S    I  CORRECTION  FOR  THE  EMPIRICAL  CALC  BETUEEN  KRAMERS 

->nn  I 

_.  w  w  ! 

213  '    PLOT    SETUP  r 

IMOIIT       "TMDIIT       1       CnO      OI    PTTCO      PO      1      CPO      rOT  "       O 
ii4iwi  A • «i    w  i        1      iwr\      ri_wili_.il      wii      c      iwii      wit  I        |  4 

TC     P-  I       TUCM     OI    PTTCO      IC      "701-7      "UDCI     " 
li         _       a        iti_,i«      ii_wiiL.li       Aw       i— 'if       i  ii     _  i_ 

T1/I  TC     n-T      TUC»I      Dl    flTTCD      TC      POT      •  T  kj  TCQUJOI    " 

*_    A  *>  _  I  _       '   l_  i     ,    ,  l_   i   i         l_.WII-.lv  _W         Ultl      |  Lll    I    LHili    iL 

"I1C       I       DC      CIIQC      Dl    OTTCD     niD      CMTTPUCC      &DC      OQPpCPl    V      CCT      .  ...         CI.ITTPU      1       1      1      T  HI     DPCTT       1 

A.    A    W  WW  WWtS_.         l_.Wtl_.lt         W  A  I  W«t  A     I    Wl   ■-.  W         i    ■  1    '  I —  lltwV_.ri_.l  wC    I  ACT*  W««   A     I    Wll  1     f    A.    •   W  Alt        "J.!     I  I 

">  ">  n»      I      CTHJTT 

A-  A.  W  W   A  14  A     I 

~>?2    '    GRAPHIC0    ON 

243  LINE    TYPE    1 

->cn>  iitcypnOT     ic     1  ~>n>     ic     on 

A-W-l  V    A   _.  *»l       Wl  I     I  AW1AA.W,AWtWV 

ten  UTMnny    n     1  n    n>       c 

-i-7«  avee     t         i     p    m    c 

— .     '     W  IWIA-W  AftAfWfWfW 

-103  pi    to    nf-C 

-ion*  pc  t  tc     i        c 

tiki  i  noc    c 

WW  XJ  l_  Wl  i  W         W 

313  '  LABEL  HORIZONTAL  AXIS 

-7-int  coo  M-f»  to  i  n» 

w  I.  W  i       wl  i         II     "V  I     W  A   0 

no  uni  ic     M     _     nn 

"i _fl  I   -OCI       M 

w-r «  u.i  i  w  _._.      ii 

T  C  OI  kj  C  V  T      M 

2E3  MOVE    5, -.37 

273  LABEL    "VELOCITY    !m/_!" 

-TOflj  I       I    ADCI        HCDTTPAI        _VtC 

WWU  _>ri  _/_.—.         WWIIIAWIIW        ll/tiw 

7on  i  npc    o 

WWW  W   W  I   V  £  w 

'DA  CPO      M_«      TP  C      CTCO  1 

41  r%  mpi  ic     —     <•     m 

AV  i   i  W  V  C  •-,(! 

423  LABEL    M 

/(TOI  kICVT       M 

-r  w  %J  i1_n   i        ii 

^^m  I    n  TD      O  T    /-> 

4cm  1  npc    c 

<CD  MPI  IC       _  1        C  ">C 

-»W  W  1    .W   I,  w.  A     >    W    I     .    A.  W 

473  LA2EL    "h    [c_»l    _-"-2    3ec"-l    C"-l]" 

<  OO  I        T  T  Tl     C       DIPT 

a  on  1   p  t  d     a 

egg  L.0RS    4 
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E 1 2 

Clffl 

533 

C  /in 
CJ  w  u 

C"7n* 

coo 

w>  \j  xj 

can 

-j  ~»  \i 

E33 
£12 

E/i  m 
TV 

E52 
ccn 

(_»  u  u 

CTf* 

COJ 

can 

"7fltPI 
■  \j  1/ 

-7  1  n 

7Z3 

im 

753 
-ico 
n» 
ion 
->on\ 
oma 
o  1  n 


LftBEL    "COMPftRISION   OF   CAMPBELL    TO   ftFSL 

urn  ipc        ± 

C£ I ZE   3 

1   AOCi       "UTTLi     pnppcnj  T  OKJ  » 

rcnc     1 
1 
u-  1  -1  —  n 

1 1_  1  -.  -  .1 

V    W    A    W    ~«/ 

Mnnc     m     n> 

1    Tiic     TYPC     c 

FOR    '»'~1    TO    12 

Do.-D-.  .  j(i  1,  /  n  /  1  n»m      wm.. 
tc     <3-^^  tnnnn    tuc>»j 


1/  1 


oma 


►.j  _    o  m  c 

GOTO   C>lc; 
END    IF 

tc     a^\*nnn     tucm. 

K  I*.  -    '74 
N  -  .  E  1 S 
GnT0    Cs 1 c 

ciin     tc 

^_l^W  A  < 

rp    p,\«j>    T  H  c  M 

K  !t  -  .  E  1 5 
m  -    4.CC 
cnTn    r  -  1  - 

WW   I    w       w  \J  a  v* 

cwn    ;c 

jc    P2.>4    THEN 


W    W  tj 

o  1  a 

W  TO 

252 
°E2 
07m 

922 

OOffl 
D  1  m 
a~>n 
a~in 


N -.  395 

CMH      TC 

1 —  '  *  ■ —       A  I 

1/ 1,  -     OOI 

N~    ,TI 


,3  1 


u  —  /  f  i,  *  1/  '  n  v^/ii^/n/inini 
MOI  IC     II-  1  -j     u-  1  -1 
i-iDai  1    1 1    u 
1 1  - 1  J— 1 1 

Hcld-M 

►JCVT      II 

I     THJC       JYPC       C 

1 1-  i  j-n 

Mnnc     n    n 

LJ  -  1    -I  _  .-» 
t  iw  a  w  — iu 

COD      I  I-  1       Tn       1  01 
r  ui »     <*  -  1      1  w      a  o 

R--c- .«!  1,  (  n  /  t  nn  wm.. 
tc    R3-<cg    THEN 

N'  'd -.  91  *Pr "  .  31  *Ra " 
ELSE 

Nt4d-.S*Pr"  . 31 *Ra"  .E 

END    IF 
\  n~in     1        u-»i,.j,r  —  ,._-*, i'/n 

l  ts    t   \j  ilt«CAwww'     1     Cwb-*f»/w 

1  non  lj-m — ui/  /n 

11mm  Mnnc     11-IJ    u-1-1 

1  t  1  n  n  p  Ck  y    11    u 


942 
acn 


oon 

w  w  t» 

992 
1  nnn 

a  w  w 

1  n  1  in 

1 223 

1  nin 

I  V  w  V 

iUTW 

1  (ICfl 

A   1/  WU 

men 


1  f»c 

A  u  w 
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i  i  tr* 

i  i  _  w 

i  i  a  n 

I  1  cm 
i  icn 
i  i  no* 

I I  S3 

*  i  -i  *j 
i  inn 

1   ">  1  f» 

*  i-  U  \J 

1  17(1 

n«ffl 

A  *_  -r  u 
t   1Cf» 

i  -icm 

1.  i-  o  t> 

i  ion 

A  i_  W  W 
1     ^f*/7l 

nm 

t  ITS 
\  IT* 

i  t  a  n 
1 252 

i  ten 

i  -on 

1  ion 

i  /(  no 

1  <  en 


I  I  _  I  A  -I  I 

Hold^H 
NEXT   V 
I 

C5I2E    3 

I     IMC      TVOC      C 

L.  A  l\C  I      I    l      4-        -I 

Mfinc     i        ccr 

i  .V/  VI-         *     ,     •    ~l*t 

npau    ->       cc 
Mnnc    i       cc 

I    I  W   V   b-         «-/    |     «    bJ  ~J 

LINE    TYPE    1 

I    ADC!         "      OflMODCI    I    It'OICTUl" 

LtlUi_L.  Ul>!   II      ULUUt    IllliL     II    I     I 

i_  imi_      i    i  r  i_     i_i 
MOIIC      C  EC 

HWVU.         W|      •    ~J  -J 

nq&u    -7       cc 
Mni  ic    o       cc 

LINE    TYPE    1 

LABEL    "    AFGL ' KRAMERS  ) " 
Mm  ic    £»   a 

I     TMC       TVOC        4 

L_    X    I   «l_  •        I      i       l_  -T 

coo    n— 1     m    c 

i     L/i  I       V   — X         I   u       o 

U-O      A  0~<CCC--Z  I  I       700T7C-/   H  I.  I    \*      B7U1C01   \ 

npai.i    ii    u 

"ill   l»I         V    ,   i    i 

N  c  X  T  ' ' 

moiic  "3 

noAii  a 

mhi  ic  c 

I     TMC       TVOC        I 

LABEL    "EXPERIMENTAL' 

c»in 
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APPENDIX  D 

COMPARISON  TEST  DATA 


isa 


i  oo 


3   Sep    1388 

~W  -Jw-fi  j 


;;■  >-•■  1 i» ■:.:;.:♦.  *j*  )&*: 


MPS 

3§  - 


14!  15 


!50 


100 


50 


TIME   RVERFIGED  CT-2 


IB' 

10* 

~    18- 

►- 
u 

10- 
10- 
10- 


10 
M: 


9G 


LOG    OF    CT-2 
3    Sep    1988  NPS 


iH#^ 


14:83  14:10  M:15 

Ti  m 


71 


150 


y  100 


50 


0 


15:00 


3   Sep    1388 



lists  ST.' 

•'V    -    ,HHH* 


NPS 


15:05  15:10 

TIME    < LOCAL) 


15:  15 


150 


100 


50 


TIME   MVERRGED  CT-2 


i  m-»  '.m-t 

LOG    OF    CT-2 

3    Sop     1988  NPS 


10' 

10' 

?    10- 
<j 

10- 

10- 

18' 


10" 


15:00 


15:05 


15:10 


Tii 


15:  15 


72 


150 
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